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ABSTRACT 


Title of Thes.is: A Study of Low Energy Cosmic Rays at 1 A. U. 

James H. Kinsey, Doctor of Philosophy, 1970 
Thesis directed by: Frank B. McDonald 

Professor of Physics, Part Time 

The results from the two scintillator AE versus E - AE telescopes 
on IMP- III and IMP- IV and the solid state telescope on IMP- IV are ana- 
lyzed and the resulting proton and alpha particle fluxes presented . 

The low energy flux time histories and energy spectra are shown for 
the energy interval 18.7 to 81.7 MeV/nucleon from June 1965 to April 
1967, and in the interval from 5.2 to 81.7 MeV/nucleon from May 1967 
to August 1968. 

A comparison of the quiet time spectra of both protons and alpha 
particles is made. It is shown that the results after solar minimum 
in 1965 do not agree with currently accepted theory in the low energy 
region of the spectrum considered. Further, it is shown that the 
reason for this may be because of a hysteresis in the particle fluxes 
with respect to energy. 

In September of 1966 a very abrupt decrease in proton fluxes at 
energies of 70 MeV and below is shown to have occurred which did not 
recover again to its previous level. This decrease is attributed to 
a change in the properties of the propagation medium following the 
large solar proton event of 2 September 1967. 

It is found that there is a fairly flat ratio of He^ to He^ + He^ 
in the energy range considered, with a value of about 77®. 



Further evidence is presented for the existence of both recurrence 
events with 27 day periods which are related to large calcium plage 
regions on the sun and co-rotating regions which produce discrete pro- 
ton events observed at earth. These observations serve further to 
establish the source of protons with MeV energies. 

A two component model of low energy cosmic rays is investigated 
in relation to the IMP-IV proton and alpha particle data. This model 
treats the observed cosmic ray flux as the sum of .two power laws in 
energy, one with a negative exponent which is taken as the solar com- 
ponent and the other with a positive exponent which is identified with 
the modulated galactic primary component. It is shown that this model 
fits the observed spectra extremely well in the energy and time inter- 


vals considered. 



A STUDY OF LOW ENERGY COSMIC 


RAYS AT 1 A, U.* 


b Y ■ 

James H. Kinsey 


Dissertation submitted to the Faculty of the Graduate School 
of the University of Maryland in partial fulfillment 
of the requirements for the degree of 
Doctor of Philosophy 
1970 


‘^Supported in part by NASA Grant NsG-695 



ACKNOWLEDGEMENTS 


It is not surprising in retrospect to find that a number of people 
have contributed to a greater, or lesser, extent in making this work 
successful. Above all, I am sincerely grateful and indebted to my 
thesis advisor, and mentor on matters cosmical. Professor Frank B. 
McDonald for letting me have -the privilege of working with him. 

Further, I am most thankful for his tolerance with my occasional • s ide 
excursions in the process of this research, which he always very tact- 
fully and philosophically directed back to the main stream. 

Drs . B. J. Teegarden and G. M. Simnett receive very special thanks 
for making available their work on detector calibration and gain shift 
measurements for these experiments, as well as for the many suggestions 
they put forth during the process of the research. Dr. D. E. Hagge is 
most cordially thanked for sharing the data output from the low energy 
solid state particle telescope, which he designed. 

Drs. V. K. Balasubrahmanyan, E. C. Roelof, R. Ramaty, D. Reames, 

T. L. Cline, C. E. Fichtel, P. A. Abraham, E. Boldt, P. J. Serlemitsos, 
and U. D. Desai all took a considerable interest in this research and 
were always willing to make -their ideas available. Their suggestions 
and criticisms were invaluable. 

Messrs. Ciro A. Cancro and Paul J. Janniche are most gratefully 
acknowledged for their design and construction of the flight electronics 
for the experiments . 

Most of the thanks for the unglamorous task of initial data pro- 
cessing must be given to the programming efforts of Messrs. P. J. Barry 
and 'C. Cattell and Mrs. J. B. Brooks. I am also most grateful to 

ii 



Mrs. G. M. Pearson for keeping track of data tapes, and handling 
miscellaneous computer housekeeping chores. 

I also wish to thank Messrs. A. Flanick, W. Davis, and C. Glasser 
for their efforts in the laboratory in preparing the detector packages 
and providing me with calibration data. 

Messrs. A. Thompson, H. Trexel and F. Shaeffer have my thanks for 
their drafting efforts in my behalf. 

My grateful thanks go to Miss N. Ahlfeldt for a very neat and 
professional job of typing. 

I want to thank Dr. D. Klinglesmith for providing me with the 
least squares fitting routine, used . The use of this subroutine called 
"Stepit," copyrighted by J. P. Chandler in 1965 and distributed by 
the Quantum Chemistry Program Exchange, Indiana University, Bloomington, 
is greatfully acknowledged. 

Finally, I cannot thank enough my wife, Nina, who has foregone 
much pleasure and convenience, provided me with invaluable moral 
support, and kept me fed and clothed that I might complete this work. 


iii 



TABLE OF CONTENTS 


Chapter Page 

ACKNOWLEDGEMENTS ii 

LIST OF TABLES vi 

LIST OF FIGURES vii 

I. INTRODUCTION 1 

A. General I 

B. Galactic Cosmic Rays 6 

C. Cosmic Ray Modulation 10 

D. Solar Proton Events 12 

E. The Scope of This Study • 13 

II. THE EXPERIMENTS 18 

A. Theory of the Detection Technique 18 

B. Csl Scintillator Telescope for IMP-III 

and IMP- IV 23 

C. Silicon Solid State Telescope for IMP-IV 25 

D. The IMP-III Experiment 27- 

E. The IMP-IV Experiment 29 

III. DATA ANALYSIS 31 

A. General Format of Data '. 31 

B. Calibration of AE Versus E - AE Response ..... 32 

C. Matrix Histogram Analysis 36 

D. Computation of Particle Flux Time Histories 

With Fine Time Resolution 39 

IV. THE OBSERVATIONS 45 

A. Time Histories of Low Energy Proton Fluxes 

From June 1965 to August 1968 45 

iv 



Chapter Page 

B. Proton and Alpha Particle Spectra and 

Cosmic Ray Modulation 50 

C. Comparison of Quietest Time Energy Spectra 

From IMP-IV with Solar Minimum Spectrum ..... 54 

V. A TWO COMPONENT MODEL FOR LOW ENERGY COSMIC RAYS 

AND ITS IMPLICATIONS 56 

A. The Model 56 

B. Least Squares Fitting of Model to Data 59 

C. Interpretation of Results 63 

D. The Relation of the Two Component Model 

To Other Observations 66 

APPENDIX A GEOMETRY FACTOR FOR A TELESCOPE 

WITH CYLINDRICAL SYMMETRY 71 

APPENDIX B COMPUTATION OF AVERAGE PARTICLE PATH LENGTH 

IN CYLINDRICAL GEOMETRY 73 

APPENDIX C THE BINOMIAL AND POISSON DISTRIBUTIONS 75 

APPENDIX D COMPUTATION OF COUNTING RATES 

BY THE DEAD TIME APPROXIMATION 82 

REFERENCES 85 



LIST OF TABLES 


Table Page 

I. Satelli-te Data . 89 

II. Coefficients .of Telescope Response Functions .... 90 

III. Analysis Parameters 91 

IV. Semi-Annual Count of 29 MeV Proton Flux Increases . . 92 



LIST OF FIGURES 


Figure 


Page 


1. Primary differential energy/nuc leon spectra of cosmic ray 
protons and helium nuclei observed near earth near solar 
minimum in 1965. Solid and open symbols represent 
measurements of protons and helium nuclei, respectively. 

1. IMP-3, June-Dee. 1965, Fan et al. (1965a) 

2. IMP-3 and balloon, June 1965, Balasubrahmanyan, 

Hagge, et al. (1966a, b) 

3. balloon, June 1964, Waddington and Freier (1966) 

4. balloon, June 1963, Freier and Waddington (1965) 

5. balloon, June, 1964, 1965, Ormes and Webber (1966) 

6. balloon, 1954, McDonald (1958) 

7. - balloon. May 1965, Hofmann and Winckler (1966) 

8. balloon, April 1963, Anand et al. (1966) 

Only data representative of solar minimum have been used. 


(from Gloeckler and Jokipii, 1967) 93 

2. Comparison of Deep River neutron counting rate with 
smoothed sunspot number for the current cycle. Mean 
of previous- 19 cycles superimposed 94 


3. Schematic of AE versus E - AE particle telescope. A 
denotes detector measuring AE and B the detector meas- 
uring E - AE. C is the anti-coincidence detector. ... 95 


4. AE versus E - AE response curves for hydrogen and helium 
isotopes, for the scintillator telescope used on IMP-III 


and IMP- IV 96 

5. AE versus E - AE response curves for hydrogen and helium 

for the solid state particle telescope on IMP-IV 97 

6. Schematic diagram of geometry of a generalized circular 
cylindrical telescope of height h and base radii of R-^ 

and R 2 *. . 98 

7. Schematic of IMP scintillator telescope assembly 99 

8. Schematic of IMP scintillator telescope showing, critical 

dimensions and composition of layers 100 

9. Geometry factor G versus penetration of particle into B 

detector for IMP scintillator telescope . . . -. 101 

10. Schematic of IMP-IV solid state detector telescope 

assembly 102 


Vll 



Figure Page 

11. Schematic of IMP-IV solid state detector telescope show- 
ing critical dimensions and composition of layers .... 103 

12. Geometry factor G versus penetration of particle into 

B detector for IMP-IV solid state detector telescope. . . 104 

13. Block diagram of IMP-III cosmic ray experiment elec- 
tronics 105 

14. Block diagram of IMP-IV cosmic ray experiment elec- 
tronics 106 

15. IMP-III AE versus E - AE matrix showing proton line and 

six energy bins: (1) 18.7-29.2, (2) 29.2-39.7, 

(3) 39-.7-50.2, (4) 50.2-60.6, (5) 60.6-71.1, and 

(6) 71.1-81.7 MeV 107 

16. Proton AE versus E - AE histograms. Region under the 

dashed curve represents the background component. The 
shaded area represents the true proton component 108 

17. (a) 24 hour average proton fluxes 109 

(b) 24 hour average proton fluxes 110 

(c) 24 hour average proton fluxes Ill 

(d) 24 hour average proton fluxes 112 

(e) 24 hour average proton fluxes 113 

(f) 24 hour average proton fluxes 114 

(g) 24 hour average proton fluxes 115 

(h) 24 hour average proton fluxes 116 

(i) 24 hour average proton fluxes 117 

(j) 24 hour average proton fluxes 118 

18. IMP-III and IMP-IV monthly averages of proton fluxes and 

monthly averages of the Deep River hourly neutron rate . . 119 

19. IMP-IV 96 hour average proton fluxes 120 

20. (a) Quiet time proton energy spectra 121 

(b) Quiet time proton energy spectra 122 

21. (a) Quiet time alpha particle energy spectra 123 

(b) Quiet time alpha particle energy spectra 124 

22. Modulation of proton energy spectra since solar minimum 

with respect to the solar minimum spectrum 125 

23. Modulation of the alpha particle energy spectra since 

solar minimum with respect to the solar minimum spectrum. 126 

24. Time history of the 24 hour relative particle fluxes 

from the IMP-IV discriminator level 6 127 


vm 



Page 


Figure 

25. Representative 96 hour proton and alpha particle energy 

spectra with the least squares fit to the two component 
function 128 

26. (a) Representative alpha particle-to-proton ratios of the 

96 hour spectra 129 

(b) Averages of alpha particle-to-proton ratios at 

discrete energies from 24 May 1967 to 20 August 1968 

showing standard deviations 130 

27. 76..4 MeV versus 5.2 MeV proton fluxes 131 

28. 5.2 MeV proton flux versus the energy at which minimum 

flux occurs 132 

29. 76.4 MeV proton flux versus the energy at which minimum 

flux occurs 133 

30. Minimum proton flux versus the energy at minimum flux . . 134 

31. The computed power law indices versus the 5.2 MeV proton 

flux 135 

32. The computed power law indices versus the 76.4 MeV proton 

flux 136 

33. 74.2 MeV/nucleon versus 5.9 MeV/nucleon alpha particle 

fluxes — — 7- .... 137 

34. ' 5.9 MeV/nucleon alpha particle flux versus the energy/ 

nucleon at which the minimum flux occurs 138 

35. 74.2 MeV/nucleon alpha particle flux versus the energy/ 

nucleon at which the minimum flux occurs 139 

36. The minimum alpha particle flux versus the energy / 

nucleon at minimum flux 140 

37. The computed power law indices versus the 5.9 MeV/ 

nucleon alpha particle fluxes .............. 141 

38. The computed power law indices versus the 74.2 MeV/ 

nucleon alpha particle fluxes 142 

39. The 24 hour average proton fluxes at 5.2 MeV for IMP-IV. 

The times large confirmed solar flares occurred are marked 
F. Central meridian passage of calcium plage regions are 
marked by the plage region number. Recurring regions are 
connected by horizontal lines 143 


IX 



Figure Page 

40. 107» variable time averages of proton fluxes for 3 months 

of IMP- IV data showing solar proton increases caused by 
co-rotating regions ... 144 

41. .6 hour integral proton fluxes for energies > 20 MeV 

from IMP- III - 145 


x 



I. INTRODUCTION 


A. General 

Until 1958, with the launching of the first scientific satellite, 
all measurements of the energetic particle population in the solar sys- 
tem came from ground based, balloon borne, or rocket borne particle 
detectors. Balloon experiments measured- the fluxes of primary cosmic 
ray particles in the range of kinetic energies from ~ 10 to 10 ev/ 
nucleon. Above this energy it was necessary to rely on the measurements 
of the secondary cascade -particles created in the earth's atmosphere by 
high energy primaries with kinetic energies of between 10^ and 10^ ev/ 

g 

nucleon using ground based air shower experiments. Beloxtf about 10 ev/ 
nucleon and also in the range between the balloon and air shower measure- 
ments, little was known of the particle populations. Artificial satellites 
quickly became extremely attractive as particle experiment platforms since 
it is possible to get completely above the last vestiges of atmosphere for 
long periods and in some cases to get outside the earth's magnetosphere. 

This is a vast improvement over balloons, which can only reach altitudes 

-2 

where there is still a few g cm of air remaining above them and only 
for periods of tens of hours. 

The motivation and justification for putting particle detection 
experiments aboard satellites is several-fold. First of all it is 
desirable to study the intrinsic particle populations at 1 A. U. free 
of as much of the effects of the atmosphere and magnetosphere as possible. 
Further, one would like to ascertain the interstellar spectrum of parti- 
cles and in order to do this it is advantageous to start with as few 
extraneous effects as possible. Finally, it is of interest to study 
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the propagation of particles coming from the sun as well as the sun as 
a source of energetic particles. 

The actual form of the interstellar particle energy spectrum is not 
directly measurable near the earth because of the solar modulation of 
particles passing into the solar system. It is also now known that the 
sun itself accelerates copious quantities of particles in conjunction 
with solar flares and at low energies is producing a quasi-continuous 
component of energetic particles. This is contrary to the previously 
made assumption that during solar quiet times the observed energy spec- 
trum of cosmic rays was predominantly the solar modulated galactic spec- 
trum. Such energy spectra for both protons and alpha particles are shown 
in Figure 1 for the period in 1965 near the minimum of solar modulation. 

A local minimum in the neighborhood of 50 MeV/nucleon complicates the' 
behavior of this spectrum in the energy interval below about 300 MeV/ 
nucleon. The problems of the form of the solar modulation and the 
amount of solar particle contribution to the lower energy portion of 
the cosmic ray spectrum must be solved in order that the’ true interstel- 
lar spectrum may be determined, as well as contributing to the understand- 
ing of the sun as a source of energetic particles and the understanding 
of the large scale structure of the interplanetary medium. 

The importance of knowing the true interstellar particle energy spec- 
tra, particularly for protons, with energies below 300 MeV/nucleon is 
quite apparent. Because of the steepness of the spectrum at higher energies 
it is possible that most of the energy transported by cosmic rays in 
interstellar space is resident in protons having kinetic energies below 
1 GeV. Cosmic rays have an energy density in space that is comparable 
to the energy- density resident in all forms of electromagnetic radiation. 
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The amount and form of this energy is most important for the equilibrium 
balance of energy in the galaxy, for instance. Balasubrahmanyan et al. 
(1967) have shown that the observed heating of interstellar HI clouds 
using an assumed rate for the heating can be accomplished by the coulomb 
interactions of cosmic rays with the matter in these clouds if the low 
energy fluxes are high enough. Theory shows that the lower energy com- 
ponent of the cosmic rays would do most of the heating. 

Since the sun plays so important a role in controlling the particle 
fluxes seen at earth, it is of interest to compare cosmic rays and 
solar activity with respect to time. Figure 2 shows the average of the 
smoothed sunspot numbers for cycles 1 through 19 superimposed on the 
current cycle number 20. On this plot is also shown the monthly averages 
of the Deep River hourly neutron rates with an inverted rate scale. There 
appears to be a time lag of about one year of the neutron rates behind the 
sunspot numbers. It has been shown by several workers that this apparent 
hysteresis between solar activity and neutron rates is probably not sig- 
nificant in itself. Simpson and Wang (1967) and Hatton et al. (1968) 
have shown that the coronal FeXXVI green line at A5303 correlated very 
well with the cosmic rays and with no apparent time lag. Balasubrahmanyan 
(1969) showed that the magnetic index Ap, which is closely correlated with 
the solar wind velocity, also correlates well with the neutron rates dur- 
ing solar cycle 19. 

Since the observations of Lange and Forbush (1942) that there was 
a correlation between sea level ion chamber counting rates and solar 
flares, it has been known that the sun was a source of energetic parti- 
cles. Although these solar particles are sometimes referred to as 
"solar cosmic rays," the term "cosmic rays" is usually reserved for 
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particles originating outside of the solar system. There is substan- 
tial evidence that non-flare associated particles are released by the 
sun in so-called recurrence events with a 27 day period as well as in 
connection with isolated solar active regions of less than 27 days 
lifetime (Fichtel and McDonald, 1967). The 27 day period is the 
equatorial rotation period of the sun. The proton flux increases 
associated with this rotation occur approximately when a large persis- 
tent calcium plage region passes central meridian each rotation. 

Briefly, then, the observed characteristics at 1 A. U. of these 
four different classes of particles are 

1. Galactic cosmic rays which make up the bulk of particles de- 
tected at energies of about 50 MeV/nucleon and above during solar 
quiet times. 

2. Solar flare particles which are a transient component with 
increases of as many as 6 orders of magnitude above quiescent 
rates at low energies. These particles show very rapid rise 
times in their fluxes followed by exponential decays of the 
order of days. The energy spectra of flare particles is steep 
and of negative slope. 

3. Recurrence event particles whose fluxes rise slowly to a 
plateau for several days and then decay again so that the flux- 
time profile is quite symmetrical. These particles are recognized 
by their 27 day recurrence pattern usually coinciding with the 
time of meridian transit of large persistent calcium .plage 

1 regions on the solar disk. 

4. Increases in low energy fluxes during quiet times similar in 
character to the 27 day recurrence events , but which occur only 
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occasionally. These particles are associated with large active 
regions near the central solar meridian. 

The interplanetary magnetic field is composed of the various 
magnetic flux loops originating in the sun and transported radially 
outward by the solar wind. Near the sun, as shown by Parker (1963), 
the mean field has the configuration of an archimedian spiral, which 
much of the time dominates the propagation of low energy particles 
moving from the sun outward. On the other hand, the irregularities 
that exist in these fields cause scattering of particles and give rise 
to diffusion- like flux characteristics. The cosmic rays coming into 
the solar system interact with these moving fields and their irregu- 
larities, and it is the form of these interactions that determines the 
form of the flux modulations. 

The earth's field' also dominates the motion of charged particles 
entering the magnetosphere. The interactions are such that at a given 
latitude there exists a critical momentum below which particles are 
deflected back into interplanetary space rather than being able to 
penetrate near to the earth's surface. This critical momentum, or 
rigidity, is inversely proportional to geomagnetic latitude. This 
effect is used in balloon latitude surveys such as carried out by 
Neher (1967) as a naturally existing magnetic analyzer to determine 
the spectrum of cosmic rays. 

Most cosmic ray experiments detect charged particles by measur- 
ing the signal produced when they lose energy in passing through the 
material of the detector. Knowing the effective area and solid angle 
subtended by the detector for incident particles and the number of 
particles counted during a specific time period, it is possible to. 
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express the observation as a specific flux, or intensity, of particles 
with units of number per unit time, unit area and unit solid angle. Such 
fluxes may be studied in terms of essentially four experimental observables 

1. Energy dependence. 

2. Charge and mass composition. 

3. Time dependence. 

4. Directional dependence 

Over the past few years there has been a several fold improvement in 
resolution of these observables because of continually improving 
experimental techniques . 

B. Galactic Cosmic Rays 

During periods of minimal solar activity when there are few, if 
any, major flares, the residual particle flux was generally accepted to 
be coming from outside the solar system. The modulation effect becomes 
small at higher energies than ~ 12 GeV and the differential particle 
flux, i.e., the flux per unit energy interval, has been found to fit 
a power law in kinetic energy per nucleon quite well at higher ener- 
gies and is usually expressed as 


dJ/dE = KE“ 7 


( 1 ), 


where the nominal value of the power law index is y ^ 2.5. The integral 
flux J(> E q ) , which is the quantity measured when using the geomagnetic 
latitude effect, is defined as the total flux above some energy E Q and 
is related to the differential flux by 


J(> E 0 ) = J* (dJ/dE) dE 

Eo 


K'E„ 


-(7-1) 


( 2 ). 
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The outstanding features of the galactic component at 1 A. U. can. be 
characterized by the following general observations: 

1. A very long term time variation with a period of 11 years, which 
is the period of the solar cycle itself. The cycle of the particle 
fluxes lags the sunspot cycle by about one year. 

2. Nearly an isotropic distribution of particle arrival directions 
is observed. 

3. A monotonic ally decreasing energy spectrum above 300 MeV which 
approaches quite well a power law in energy with an index of -2.5 
above about 12 GeV. 

4. A composition which differs quite markedly from universal 
abundances 'in that there are more heavy nuclei than would be 
expected . 

-3 

5. The energy density of cosmic rays is approximately .5 eV cm 
which is comparable to other forms of energy in the galaxy such 
as that contained in electromagnetic radiation and the inter- 
stellar magnetic fields. 

6 . Frequent decreases in counting rates of a temporary nature 
associated with magnetic storms and frequently following solar 
flares and closely coinciding with the arrival of the resultant 
plasma wave at the detector. These decreases are called Forbush 
decreases . 

Any model for the production of cosmic rays must invoke sources 
energetic enough to account for the total energy observed and must also 
provide a means of transforming such energy into kinetic energy of the 
particles. Several source mechanisms have been considered in the lit- 
erature. One of the earliest such mechanisms discussed was the Fermi 
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acceleration mechanism, a process in which particles could increase 
their energies by multiple reflections from moving magnetic fields. 

A very clear discussion of this mechanism is given by Morrison (1961). 
One of the most attractive sources up to the present has been the super- 
nova, which provides adequate energy and has a frequency and spatial 
distribution in the galaxy, such as to account for the isotropy observed 
in cosmic ray fluxes, as well as being a potentially good source of 
heavy nuclei which are abundant in the cosmic ray primaries. Supernova 
acceleration models have been discussed by Burbidge, Burbidge, Fowler, 
and Hoyle (1957), Colgate and Johnson (1960), Colgate and White (1963), 
and Arnett (1966). Most recently since the discovery of pulsars there 
has been some discussion that these objects could be the source of 
cosmic rays. Of the above source mechanisms, the Fermi process is 
probably the least likely since the energy spectrum of particles that 
can be calculated using this model is steeper by at least an order of 
magnitude than is observed. 

The nuclear composition of cosmic rays and their comparison to 
universal abundances is a very important source of information about 
not only the source of the particles themselves, but also about the 
propagation of these particles through the interstellar medium and the 
nature of this medium itself. The abundances of heavy nuclei with 
respect to protons is much higher in cosmic rays than in the universal 
scale. This fact can be interpreted either as being due to sources 
rich in heavy elements, which would be the case for supernovae since 
they have evolved their element burning up to iron, or the superabundance 
of heavies can be attributed to preferential acceleration of heavies 
over the lighter nuclei (Ginzburg and Syrovatskii, 1964). The latter 
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is not likely on the basis of solar particles which are deficient in 

heavies. Further, there is an anomolous abundance of the light elements 

He , Li, Be, and B in cosmic rays. These are most easily accounted for 

by fragmentation of the heavy nuclei on the interstellar hydrogen in 

their propagation from the source through the interstellar medium to 

the observer. Assuming that the source abundances of all components 

are the same as the universal abundances, it can be shown that the 

superabundances of this light component can be accounted for if the 

integrated hydrogen along the path of propagation is of the order of 
-2 

4.1 g cm (Ramaty and Lingenfelter, 1969). 

Very closely related to .the nuclear content of cosmic rays is the 
electron component and its relation to the electromagnetic radiation 
associated with cosmic rays. One might expect to find equal numbers 
of electrons and positive ions present in cosmic rays. On the contrary 
the ratio of electrons to protons observed in cosmic rays is about 1%. 

For a time it was believed that the electrons might be produced locally 
in interstellar space by the high energy interaction of cosmic rays, 
such as p - p reactions and neutron decay (Ramaty and Lingenfelter, 

1966). L* Heureux (1967) has shown that this theory is not consistent 
with observed results, because the spectra of electrons and positrons 
produced by this mechanism fall well below the observed spectra. Further, 
the collision theory predicts that there should be an excess of positrons 
in the cosmic rays. That this is indeed not the case has been shown 
(Hartman, et al. , 1965). One is left then with the conclusion that 
electrons are probably accelerated by a hierarchy of sources. Evidence 
that this is true is exhibited in supernovae remnants such as the Crab 
Nebula, which is a strong emitter of polarized radio emission consistent 
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with synchrotron emission from electrons trapped in the remnant's magnetic 
fields. Recent advances in x-ray and gamma ray astronomy indicate an 
abundance of sources for these forms of radiation in the galaxy (Fazio, 
1967; Morrison, 1967). The assumption is that these observed sources of 
photons are probably where the cosmic rays are being accelerated, or at 
least part of them. 

C. Cosmic Ray Modulation 

Any theory that would attempt to infer properties about interstellar 
space and the galaxy in general from the observed cosmic rays in the 
vicinity of 1 A. U. must include a means to account for the changes 
in the cosmic ray fluxes. Observed cosmic ray fluxes are modulated 
■with a period of 11 years. Webber (1967) has summarized and compared 
the observations of a number of workers very well up to the minimum in 
solar modulation in 1965. The' general assumption has been that in the 
absence of solar flares and recurrence events, the quiet time observed 
particle fluxes are primarily the modulated galactic component (Gloeckler 
and Jokipii, 1966). There has been to date a small amount of evidence 
that at very low energies there might be a quiet time solar contribution 
(Meyer and Vogt, 1963; Fan et al., 1968, 1969). Neher (1967) by means 
of latitude surveys of ionization rates with balloon borne ionization 
counters has set an upper bound of approximately a factor of 4 to the 
change in the integral flux of particles for energies > 100 MeV between 
solar maximum and minimum. Vogt (1962) and Meyer and Vogt (1963) found 
primary proton spectra in 1960 and 1961 following the previous solar 
minimum which were considerably higher below 100 MeV than any of the 
other observations to date. Furthermore, they found a relative minimum 
at about 100 MeV as compared to the one at about 50 MeV. There has 
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been speculation as to whether the present spectrum could return to 
values such as they observed. 

Parker (1958, 1963, 1966) has treated modulation as both a con- 
vective process whereby the scattering centers in the magnetic fields 
are being strongly affected by solar activity, and in terms of adiabatic 
deceleration of incoming particles due to the effectively expanding 
magnetic field in the solar wind. Theoretical treatments of solar 
modulation have been done by Gleeson and Axford (1968), Jokipii (1967, 
1968), and Jokipii and Parker (1967). A good summary of the critical 
aspects of the various theories with respect to the observations is 
given by Webber (1968). 

Parker's theoretical relation for modulation (1963) may be expressed 


as 


dj e (P,t) = dj oo (P)e 


1 


J 


a. u . 


w dr 

K 


(3), 


where the differential flux at the orbit of earth di is a function of 

e 

particle rigidity and time, dj^ is the unmodulated flux beyond the 
solar system and w and (c are the solar wind velocity and diffusion 
coefficients, respectively. Gloeckler and Jokipii (1966) have shown 
that the diffusion coefficient is separable into the produce of a 
function of rigidity and velocity and a function of radial distance 
and time, thus 


sc = f(P,P) g(r,t) 


(4). 


Using this last relation it is possible then to express (3) in the form 


d J e 0?>t) 


dj ro (P) 


-Tl(t)/f(P,p) 


(5). 



Following Nagashima et al. (1966), it is convenient for comparison of 
spectra at different epochs t and t' to take the logarithm of the 
ratios of the two fluxes using equation (5) to get 
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dip (P,t) _ ATI 

dj e (P,f) f(P,0) 


with 


AT] = TKf) - Tl(t) 


( 6 ), 

(7). 


Webber (1962) shows that using the relation in equation (6) with 
t corresponding to the solar minimum in 1965 and t 1 for the fluxes 
prior to this time that the form of f(P,(3) is Pj3 for P > 1 GV. The 
• 1/P0 and l/{3 dependency of the logarithmic ratio of fluxes holds for 
both protons and alpha particles; however, the modulation constant 
for protons is larger. For the time period considered, Webber shows 
that the observations support a mixture of both the diffusion-convection 
and the energy loss theories. 

D. Solar Proton Events 

As mentioned above, there are large bursts of particles following 
solar flares. These particles come as a rapid increase in the low 
energy fluxes and occur shortly after a large solar flare reaching flux 
values at the lowest measured energies of as many as 6 orders of magni- 
tude above quiet time levels. These large increases in counting rates 
are caused primarily by the proton component and hence large particle 
producing flares and the attendant high proton fluxes are often referred 
to as "solar proton events." Such events are characterized by the 
following features: 

1. An intense brightening of hydrogen alpha emission in the 
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vicinity of a solar active region in the solar chromosphere. 

This brightening is almost always accompanied, or preceded, by 
a Type IV solar microwave burst. 

2. Production of x-rays of several kilovolts energy in the same 
region as the optical flare. 

3. Arrival of protons of relativistic energies at the earth 
within a few minutes following the optical and radio flare. 

4. A very rapid increase in counting rates for all lower energies 
until a maximum counting rate is reached some hours after the 
flare. The arrival times of the particles of different energies 
usually show a dispersion in velocity, but there are also some- 
times complex variations in the counting rates that may last for 
days. In some events the low energy particles show no increase 
in counting rate until after the arrival of the flare blast wave. 

5. Following the peak flux an exponential decay in the counting 
rates commences with pre-flare levels being reached in several 
days to a week. 

6. A monotonically decreasing rigidity spectrum which can be 
fitted by an exponential in rigidity over much of its range 
(Fichte! and McDonald, 1967). 

-3 ”1 

7. Alpha particle to proton ratios of the order 10 to 10 

with the ratio decreasing with energy (Fichtel and McDonald, 1967). 


The Scone of This Stud-y 


In the following chapters it is intended that an analysis of the 


proton and alpha particle fluxes obtained from the Goddard cosmic ray 
experiments aboard the satellites IMP-III and IMP- IV will help illuminate 
some of the problems discussed in the foregoing. Because of the nature 
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of the experiments themselves and their times of execution, the following 
constraints were imposed on the analysis: 

1. The times considered were the rising portion of a solar cycle, 

going from mid-1965 to mid-1968. 

2. The kinetic energies were in the range of 4 to 80 MeV/nucleon. 

3. The charge dependence extended only as far Z = 2. 

4. The radial dependence from the sun as origin was fixed at 1 A. U. 

and .there was no directional information available. 

In Chapter II a discussion of the experiments is given in which both 
the theory of operation and the practical aspects of these particular 
experiments are considered. The general operating parameters of the in- 
dividual detectors are indicated and the anomolies of the experiments 
during their lifetimes are discussed. 

A discussion of the techniques used in analyzing the satellite data 
is presented in Chapter III I Here also the problem of calibration of 
the detectors is considered. Most of the mathematical derivations are 
carried out in the appendices. 

In Chapter IV the results of this analysis are presented in two 
forms : 

1. time histories, and 

2. energy spectra. 

The time histories are examined for regularities and anomalies during 
the three year period immediately following the minimum of solar activ- 
ity. It is found that on the basis of a number of large particle events 
observed, the most active times was in late 1966, and early 1967, while 
the quietest periods were early 1965 and late 1967. 
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Upon examination of two of the higher energy proton fluxes from the 
scintillator telescope on a monthly average basis, an anomalously large 
decrease in flux level was noted occurring in September 1966. It was 
noted that the higher energy flux was depressed by a significant amount 
more than the lower fluxes. Also the flux after the large decrease 
never returned to its pre-flare values. Further, these same monthly 
fluxes, which were at about 50 and 71 MeV, definitely showed a lag of 
about 1 to 2 months behind the sea level neutron rates, indicating some 
amount of hysteresis of the particle fluxes with respect to energy. 

An examination of the lower energy fluxes, at 8 and 24 MeV, con- 
sidered in the experiments show that the lowest energy fluxes vary by 
a greater amount than those fluxes at higher energies. Considered with 
the foregoing results of the monthly higher energy fluxes, at 50 and 
71 MeV, a positive dependency on energy, this opposite behavior of the 
lower energy fluxes seems anomalous. As it turns out, these results are 
the basis for later consideration of the lower and higher energy fluxes 
as being of different origins. 

The long term proton and alpha particle spectra are considered. 

They are compared on the basis of the diffusion-convection formalism. 
Whereas pre-1965 spectra agreed with theory in that they showed an, 
exponential modulation with the exponent proportional to 1/3, the results 
of . the present study show a behavior that is proportional to 8 , or P 
which is equivalent in the range considered, for protons. The modulation 
of the alpha particles, though present, did not fit any simple model. 

It is suggested that a possible solution to the anomoly of the proton 
modulation would be that the power spectrum of irregularities in the 
interplanetary magnetic field might be a function of time such that 
AM(K) = 1/K 3 . 
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Chapter V constitutes the formulation and examination of a two 
component model of cosmic rays at 1 A. U. to explain the difference in 
behavior of the fluxes above and below the relative minimum observed 
in most quiet time low energy spectra that have been presented to date. 
The operationally convenient fit of the sum of two power laws in energy, 
one with positive exponent, the other with negative, is made. The 
tentative association of the positive power component -with modulated 
galactic primary cosmic rays, and the negative power with solar p'arti- 
cles is made. Statistical analysis of the results of fitting the ob- 
served fluxes to the model agree remarkably well with the predictions 
of the model. 

The picture that finally emerges from the analysis of the data 
on the basis of the two component model, is that indeed the low energy 
fluxes up to the observed relative minimum in the spectrum are pre- 
dominantly of solar origin. Above this minimum, where the level of 
solar activity is low enough that the minimum can be observed, the 
spectrum is taken to be all of galactic primary origin. The short time 
resolution of 4 days used in the analysis is very much shorter than the 
11 year modulation of cosmic rays. Because of this, what is seen is 
that the galactic component is modulated very little, while at the same 
time the solar component is changing fairly rapidly and over many orders 
of magnitude. If found to be correct through further testing as the 
solar cycle proceeds, this model provides a means at low energies of 
separating the effects of the solar particles from galactic' primaries. 

In support of the sun as a quasi-continuous source of energetic 
particles of energies of about 10 to 20 MeV/nucleon, it was also under- 
taken to show that very good correlation of non-flare related proton 



flux enhancements with solar features having both a periodic property 
and magnetic activity was possible. It is shown that a number of 
reasonably frequent events of this nature did occur, and that these 
kept the low energy particle fluxes up to significant levels. 



II. THE EXPERIMENTS 


A. Theory of the Detection Technique 
When an energetic charged particle passes through matter it inter- 
acts with the atoms of the material losing energy by ionization of these 
atoms. The theory of the' collision processes involved has been worked 
out in great detail by various workers. A summary of the literature on 
the subject is given by Rossi (1952). Using the resulting stopping 
power formula for particles much heavier than the electron (Barkas and 
Berger, 1964) one may write for the energy lost per unit path length 
by a particle of kinetic energy E and charge z. 


dE _ 2rtNZe 4 z s { An? 2 8 3 - — - 6 

dx A3 2 me 3 [ in (1-3 3 ) 2 ! 2 " Z 


( 1 ) > 


where N is Avogadro's number, Z is the atomic charge of the material, 
e the electron charge, A the molecular weight of the material, -m-the 
electron mass, and 3 is the velocity of the particle relative to c the 
velocity of light. The quantity I represents the adjusted effective 
ionization potential of the material. In theory this latter quantity 
can be derived from weighted excitation energies, but in practice it 
must be obtained by performing actual stopping power, or range, experi- 
ments. The empirically derived results of such measurements may be. 
expressed as functions of the electronic charge of the material Z such 
that 


I = 12Z + 7 ; for Z ^ 12 


(2a) 


I = 9.76Z + 58.8Z -0 ’ 19 ; for Z > 12 


(2b), 
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in which the units are eV. The remaining two terms inside the brackets 
are the shell correction term C/Z and the density correction term 6. 
These two quantities are second order effects and are also determined 
in most cases from experiment. 

From equation (1) it is readily seen that for a given material the 
rate of energy loss is a function only of the particle charge z and its 
velocity 3. It is therefore possible to write the- general relationship 

dE/dR = z 3 4 (P ) (3), 


so that once the stopping power for one species of particle in a 
material is known as a function of velocity, those of different charges 
are easily found by comparison. Equation (3) may be rewritten as 


dE/m = 
dR/m 


z 3 4 (P ) 


z 3 ^ (E/m) 


(4), 


which may be integrated to give the following relation for the range 
of the particle: 


R = 2L- $j(E/m) (5). 

z 3 

Below approximately 1 GeV/nucleon it is- possible to represent E/m as 
a power law in R so that 


E = kz 3 n m x _n R n 


( 6 ) , 


with n $ 0.6 depending on the properties of the material. In general 
then one can write the approximate relation for two particles a and b 
of different charge and mass passing through a particular material as 

•t 

E b = (zb/z a ) 3n (mb/m a ) 1_n E a 


(7), 
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by which the two particles may be compared. If particle a is taken to 
be a proton then 

E = z 2n m 1_n Ep roton (8) , 

which may be used as a scaling relation to scale range, stopping power, 
or energy loss of a proton passing through a particular material to 
other nuclei passing through the same material. 

As suggested in the last chapter the- charge, mass, and energy 
distribution, as well as the arrival directions of cosmic rays are 
quantities of very great interest. The device that has proven to be 
of the greatest utility in studying low energy cosmic rays has been 
the dE/dx versus E telescope. This is a device which measures two 
parameters for a given charged particle passing through the device 
which make it possible to determine the energy of the particle and its 
charge and mass. The telescope is so constructed that a particle that 
obeys the proper detection criteria is known to have entered within a 
narrow cone defined by the geometry of the instrument. A schematic of 
such a device as used to detect cosmic rays in the experiments reported 
on herein is shown in Figure 3. The basic features of any such tele- 
scope are the three separate detectors labeled A, B, and C in the dia- 
gram. 

Consider a particle having a trajectory such as the one labeled 
2 in the schematic. The range of this particle would be greater than 
the thickness of A, but less than the thickness of A and B taken to- 
gether. Such a trajectory may be termed a stopping trajectory as 
opposed to penetrating trajectories such as the tracks labeled 3 and 
4. Detector A is usually about 5 - 10% of the thickness of B. As a 
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charged particle of kinetic energy E penetrates detector A it deposits 
an amount of energy AE. It then enters B and comes to rest giving up 
the remainder of its kinetic energy therein of the amount E - AE. The 
energy deposited in each of these detectors produces electrical signals 
which may be amplified and analyzed. For example, if A and B are made 
of scintillator material the number of photons emitted will be propor- 
tional to the energy deposited and hence can be detected and amplified 
by means of photomultipliers looking at each detector. 

To show that the measurement of AE and E - AE can indeed determine 
the charge and mass of the particle note that equation (6) may be solved 
for R as a function of E such that 

R(E) = Kz" 2 ml _ qE (9), 

where 

q = 1/n (10), 

and K is a constant for a particular detector material. The range of 
a particle of incident energy E after losing ah amount AE in A is just 

R(E-AE) = Kz“ 2 m 1 " q (E-AE) q (11). 

It is also true then that if detector A has a thickness Ax that 

R(E) = Ax + R(E-AE) (12). 

Equations (9) and (11) may be substituted into (12) and solved for AE 
as a function of E - AE with the result 

AE = [(Axz 2 /Km 1 -’ q ) + (E-AE) q ] 1/q - (E-AE) (13). 
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It is easily seen that equation (13), for a given particle tele- 
scope, represents a unique and distinct function for each set of values 
of z and m. The plots of several of these relations are shown in 
Figures 4 and 5. The particle "lines" in Figure 4 are relatively 
flatter for small values of E - AE than might be expected for the 
idealized telescope depicted in Figure 3. The reason for this is that 
in actual practice there are other layers of absorbers present in the 
construction of such telescopes in front of all of the detectors which 
are unavoidable because of the need for mechanical strength of the 
assembly. Energy is deposited in these materials as well as in the 
detectors themselves, but this energy is not detected. The lines for 
hydrogen and helium in Figure 5 are closer to theory since there is 
not as much extraneous material through which the particles must pass 
in the telescope for which these curves were plotted. In both of these 
figures the solid curves represent the AE versus E - AE responses for 
stopping particles, while the dashed curves show the continuation of 
the response for penetrating particles. 

In order to select only those particle events that correspond to 
stopping particles .coincidence circuitry is designed such that the 
desired events may be designated by the logical requirement ABC. In 
Figure 3 it is seen that of the four tracks shown only track 2 would 
satisfy this logical criterion. By restricting the analysis of events 
to those of stopping particles the following points are noted: 

1. The energy falls within the limits defined by the thicknesses 
of the two detectors A and B. 

2. Analyzed particles must enter the telescope within an aperture 
cone defined by the diameters and separation of the detectors A and B. 
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3. Spurious events are kept to a minimum by the anticoincidence 
requirement on detector C. 

4. The particles of different z and m will be easily distinguishable 
because of their different lines. 

It should also be noted that for either stopping, or penetrating parti- 
cles continuous in-flight calibration of the detectors and their assoc- 
iated electronics is provided because of the unique "endpoint" existing 
for each stopping particle. This endpoint is defined by the materials 
and geometry of a particular telescope and must represent a unique 
incident energy for a given species of particle. 

B, Csl Scintillator Telescope for IMP-III and IMP-IV 

A mechanical layout of the Csl scintillator telescope used in both 
the IMP-III and IMP-IV experiments is shown in Figure 7. As viewed, an 
incident stopping particle would enter the telescope from the right 
penetrating the light baffle and passing into the A detector labeled 
AE scintillator in the figure. The particle would continue through 
the second light baffle and stop in the B detector- labeled as the 
E - AE scintillator. Note that each of the three photomultiplier tubes 
is coupled through the open regions adjacent to each' of the three de- 
tectors. As shown the assembly is very compactly designed with assoc- 
iated high voltage power supplies and pre-amplifiers attached to the 
photomultiplier housings. 

In order to accurately compute the response curve AE versus E - AE 
as shown in Figure 4 it is necessary to take into account all of the 
material in a particle track through the telescope up to but not includ- 
ing the G detector. A diagram showing the dimensions and composition of 
the scintillator telescope used is shown in Figure 8. The aluminum and 
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plastic foam sandwich shown preceding the A detector and the one between 
the A and B detectors are primarily responsible for the flattening of 
the curves at small energies as mentioned earlier. 

To compute the response curves it was assumed that an isotropic 
flux of particles would' enter the front end of the telescope. Such 
particle tracks would generally be inclined to the axis of the 
telescope. Ideally in order to find the average response curve the 
energy losses for a large number of tracks of varying obliquity should 
be computed and the average of these energy losses for each value of 
incident energy computed. This is however not a practical approach 
since the effect of oblique tracks is second order. Rather it is 
sufficient to compute the average slant angle for a track of an iso- 
tropic distribution of entrant particle directions and use the path 
lengths so computed through the various layers of material in the range 
energy calculations. The computation of the factor by which all of the 
telescope material thicknesses must be multiplied is outlined in Appendix 
B. For the scintillator telescope in question the average correction for 
oblique tracks was computed to be l/<cos 6> = 1.02. Using the semi- 
empirical relations for range and energy of Barkas and Berger (1964) in 
a digital computer program, the computations of the energy loss in the 
various materials indicated in Figure 8 including the detectors were 
carried out. The results are those plotted in the response curves in 
Figure 4. The -thickness values used in these computations were those 
obtained by multiplying each measured thickness by the obliquity factor 
given above. 

Using equation (11) in Appendix A the geometry factors for varying 
penetration depths of tracks into the B detector were computed. These 
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are plotted in Figure 9, which has specific incident energies marked 
along the curve indicating the incident energy of a particle having that 
geometry factor. 

The design of this telescope was such that the energy/ nucleon range 
for stopping particles would nominally be between 20 and 80 MeV/nucleon. 

The actual telescope had a threshold for particle detection in the ABC 
mode of 18.7 MeV/nucleon and the cutoff for the C detector anticoinci- 
dence occurred at 81.6 MeV/nucleon. 

C. Silicon Solid State Telescope for IMP-IV 

In order to reach lower energies than was practicable with a scin- 
tillator telescope, it was decided to include on the IMP-IV satellite a 
solid state telescope using surface barrier type silicon detectors for 
the A and B detectors in addition to the dE/dx versus E telescope described 
in the last section. A schematic of the mechanical design of this tele- 
scope assembly is shown in Figure 10. The nominal values of the kinetic 
energy per nucleon for stopping particles in this telescope were 4 to 20 
MeV/nucleon. In conjunction then with the scintillator telescope described 
in the last section, such an experiment would span the region of the . 
spectrum shown in Figure 1 through the low energy minimum. 

As indicated in Figure 10 -two solid state detectors were required 
for the B detector. The reason for this was simply the unavailability 
of a reliable silicon surface barrier detector of the required thickness 
whereas it was possible to obtain detectors of half the required thick- 
ness that would meet the specifications. These two detectors then could 
have their signals added and were in this way equivalent to a single 


detector. 
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The A and B detectors were mounted in the two ends of a cylindrical 
plastic scintillator with a plastic plug of the same material after B. 
The signals from the silicon detectors were fed into charge sensitive 
preamplifiers since the number of charge pairs created by an ionizing 
particle would be proportional to the energy it deposited. The plugged 
end of the C detector was coupled to a photomultiplier. As in the 
case of the scintillator telescope described above the ABC logic was 
invoked to discriminate against all but particles coming to rest in 
detector B after having passed through A. 

The composition and geometry of the solid state telescope are 
shown in Figure 11. Note the titanium foil light shield in front of 
the A detector because of the sensitivity of silicon detectors to 
light. Also note that each of the silicon detectors has a layer each 
of gold and aluminum on front and back surfaces, respectively. The 
average obliquity factor for this assembly was l/<cos 0> = 1.04. With 
this factor applied to the thickness of material shown in Figure 11 
the range energy program mentioned in the last section was used to 
compute the response curves in Figure 5. 

In Figure 12 the geometry factor for this telescope is plotted 
as a function of the penetration depth into the B detectors, which are 
designated Bp and B2. Note that there is a pronounced jump in the 
value of G in going from Bp to B 2 - This latter discontinuity in G is 
because of the separation in the two elements of the B detector. Along 
the curve itself are indicated several incident energies for which the 
corresponding values of G occur. 

The threshold energy for the actual telescope was 4.2 MeV/nucleon, 
while the cutoff energy was 19.1 MeV/nucleon. 
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In order to provide a means of measuring even lower energy proton 
fluxes, particularly those originating in solar flares, an 8 level 
integral analyzer was used independently from the above dE/dx versus 
E arrangement to analyze the signal from the A detector. The maximum 
energy for a stopping proton in A was 4.1 MeV incident on the front of 
the telescope. The discrimination was set such that level 3 was the 
lowest that would admit a signal, that being equivalent to about 3.32 MeV. 
Level 6, which was three steps lower than 3, was set so that the integral 
flux of particles would be in the range of 1.1 to 20 MeV. 

D. The IMP-III Experiment 

The IMP-III satellite was launched on May 29, 1965. Some of the 
vital statistics for this satellite are given in Table I. The scintil- 
lator telescope was mounted so that its axis was perpendicular to the 
spin axis of the satellite. A block diagram of the electronics associ- 
ated with the experiment is shown in Figure 13. The function of the 
circuitry of interest to the present study is that part which is related 
to the in-flight analysis of ABC type events and a sample of the count- 
ing rate of these events. 

The satellite telemetry operated in such a manner that a complete 
commutation cycle was completed every 5.461 minutes. Each cycle was 
divided into. 4 sequences of 81.92 seconds duration, of which only the 
first 3 were utilized by the cosmic ray experiment. There were 2 
accumulator gates open for 39.36 seconds in each of these sequences 
giving a total of 6 accumulators effectively per cycle. One such 
accumulator each cycle was available for recording the number of ABC 
events thus providing a sample counting rate for stopping particles. 

Also during each of the 3 sequences per cycle two gates of 34.24 seconds 
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duration were open for pulse height analysis of ABC events. When the 
first ABC event occurred after one of these gates was open two separate 
512 pulse height analyzers analyzed the AE and E - AE signals from the 
A and B detectors, respectively, and stored the resultant channel 
numbers in the 34. 24. second accumulator. During the remainder of time 
the gate was open, further pulse height analyses were inhibited, thus 
effectively limiting the number of pulse height analyzed events to a 
maximum of 6 every cycle. 

The IMP-III experiment functioned well throughout the lifetime of 
the satellite. Because of temperature changes in the satellite, there 
was a certain amount of drift in the gain of the two detectors and their 
associated electronics. This drift could be compensated for quite satis- 
factorily in the analysis of the data on the ground by noting the shift 
in the channel numbers of the proton curve endpoint and computing a 
gain factor for each of the two parameters over some specified time 
period. There were at least two occasions when the scintillator crystals 
were temporarily saturated by very high particle fluxes. This occurred 
during the proton event of September 2, 1966 and again during the event 
of January 28, 1967. IMP-III effectively died in early May 1967, al- 
though the useful data available from this experiment reported in this 
study continues only until April 10, 1967. 

E. The IMP- IV Experiment 

IMP-IV, a somewhat improved version of the IMP series, was launched 
on May 24, 1967 almost providing complete coverage with IMP-III. The 
other pertinent flight data on this satellite is shown in Table I. The 
cosmic ray experiment aboard this satellite again utilized the previously 
discussed scintillator telescope to cover the nominal energy range of 
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from 20 to 80 MeV/nucleon. In what follows this telescope will hereafter 
be referred to as the MED for Medium Energy Detector. Besides this 
experiment the solid state telescope described in Section C was oper- 
ated in conjunction with it. This solid state telescope will be desig- 
nated LED hereafter for Low Energy Detector. The LED extended the energy 
range measured down to about 4 MeV/nucleon. On IMP-IV the MED was parallel 
to the satellite spin axis in contrast to its counterpart on IMP-III. The 
LED was mounted so that it looked parallel to the equatorial plane of 
the satellite, i.e., perpendicular to the spin axis. The experiment 
flight circuitry for these experiments is illustrated in block form in 
Figure 14. One significant . difference in circuit design for this experi- 
ment was the use of 1024 channel pulse height analyzers to provide greater 
dynamic range in analyzing the ABC events, although the total number of 
channels was not available because of saturation effects. 

As may be seen from the block diagram, there were two separate 
detector packages, one each for the LED and MED. A third- package 
contained the combined flight circuitry for the two telescopes. In 
this experiment a greater degree of flexibility was designed into the 
electronics so that both the LED and MED could be operated in .-a loxxr 
and a high threshold mode. The purpose for the latter change was so 
that in the high threshold mode the AE threshold would occur some 20 
channels higher than in the low mode. This would permit more alpha 
particles to be counted during high flux rates x?hen the normally high 
counting rates in the lower channels would effectively mask these counts. 
Also the LED circuitry provided for half the pulse height analyses per 
cycle being done on AB type events, i.e., on penetrating particle 


events. 
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The timing of a cycle for IMP-IV was such that there were 32 frames 
per commutator cycle of 5.12 seconds duration each giving a total cycle 
time of 2.731 minutes. Of the 5.12 seconds frame duration for 4.48 
seconds the accumulator gate was open with the remaining 0.64 seconds 
utilized for readout. There were 8 each pulse height analyses for the 
LED low, LED high, MED low, and MED high modes. The time the pulse 
height analyzer time bins were open was the 4.48 second accumulation 
time. Since the AB mode was operated for the LED half the time, there 
were then 4 LED and 8 MED ABC pulse height analyses per commutator cycle 
in each of the two threshold modes. Greater signal to noise ratios for 
these experiments was provided by causing each accumulator to be read 
out twice before resetting in order to check consistency. One accumula- 
tion sample of ABC counts for 4.48 seconds was made, one for each thresh- 
old mode for the LED and the MED each cycle. Compared to the IMP-III 
experiment, it is seen that the MED provides 16 pulse height analyses 
during the same time 6 could be performed on the earlier experiment, 
provided the counting rate was high. 

The IMP-IV experiments performed very well for the most part. Un- 
like the previous IMP experiments, including IMP-III, very little gain 
drift was noted during the effective lifetime of the satellite. Analysis 
of the data did show that the C provision of the LED part of the experi- 
ment failed in March of 1968. This failure has been postulated to have 
been caused by a cracked photomultiplier tube following the thermal 
shock of passing through the long period of the earth's shadow, which 
occurred once a year. Otherwise the experiment appeared to be function- 
ing quite well up through the last available data acquired for this 
analysis which was up to August 20, 1968. 



Ill . DATA ANALYSIS 


A, General Format of Data 

The raw data from current satellite experiments is in a form which 
is most suitable for high rates of data transmission from satellite to 
ground stations at high signal to noise ratios. This 4 bit PFM code is 
not intended for data processing in the form received, but must be con- 
verted to standard digital codes for the computers on which analysis is 
to be carried out. In the initial phases of data processing of the 
taped data as it is received from the tracking stations, this conversion 
must be carried out as well as the process of separating these data from 
different experiments and outputting this on to separate tapes for use 
by the individual experimenters, a process called decommutation. 

Several rather tedious steps of data processing had to be performed 
on the decommutated tapes before these could be used in any analysis 
schemes. Because all of the data pertinent to an experiment that was 
received by ground stations was included, there was a certain amount of 
overlap and redundancy present in the data at this initial stage which 
had to be eliminated. The result of this process was a non-redundant 
time-ordered sequence of data. Other procedures performed on the data 
in this processing (pre-analysis) stage were the insertion of pertinent 
time and orbital information where required and the reorganization of 
the format and condensation of the data. The output tapes, in a form 
for analysis, were 800 bit per inch 9-track tapes for use on IBM series 
360 computers. 

For both IMP-III and IMP-IV the taped data for this study was avail- 
able to the experimenter in two basic formats, which can be designated 
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here as (1) final data tapes and (2) matrix tapes. The former cons-isted - 
of records corresponding to the time of one commutation cycle of the 
satellite telemetry system. Each record contained the U.T. time and 
date, satellite distance from earth, and various relative coordinates 
for the particular record. The cosmic ray experiment data in each rec- 
ord consisted of the readouts of the various experiment accumulators 
corresponding to the number of counts occurring for each logic require- 
ment, e.g. ABC for one mode of the AE versus E - AE experiments, as 
well as the readouts of the pulse height analyzed events. 

The matrix tapes were the result of a form of intermediate pro- 
cessing’ which condensed the data considerably. The matrix tape con- 
sisted of sequences of n successive records which contained the pulse 
height analyzed data in a compressed matrix form, as well as accumu- 
lator sums for a specified period, usually one complete orbit of the 
satellite. By reading any n record sequence, a complete matrix of the 
AE versus E - AE ABC events could be reconstructed for a particular mode 
of any given orbit, the number in a particular location being the number 
of events having the corresponding AE and E - AE channel number. 

B. Calibration of AE Versus E - AE Response 

As shown in Chapter IX, the expected response of the cosmic ray 
telescopes was computed using the range-energy data of Barkas and Berger 
(1964). The theoretical response curves so obtained are plotted in 
Figures 4 and 5. The collective results from several orbits of pulse 
height analyzed data from one of the cosmic ray experiments can be 
interpreted in terms of one of these response curves in the following 
manner. This data can be arranged in a two dimensional array where the 
x and y axes represent the E - AE and AE channel numbers, respectively. 
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so that the number of counts occurring for a particular value of E - AE 
and AE is the value of the x,yth element of the array. Because of the 
various statistical fluctuations in energy losses of the detected parti- 
cles, there will be a maximum distribution of counts in a region of the 
array corresponding to .a particle line. 

Figure 15 shows a computer printout of such an array, or matrix, 
for the period from 4 June 1965 to 26 October 1965 for IMP-III. The 
central curved line is drawn through the maximum of the distribution 
corresponding to the proton line in Figure 4. In order to establish 
a correspondence between the physical channel numbers in the x and y 
directions and the incident energy of a proton, it was first necessary 
to know the corresponding relationship between a detector's pulse 
height analyzed channel number and the kinetic energy deposited in 
that detector. 

Scintillators exhibit a non-linearity in light output versus energy 
deposited at energies of the order of 10 MeV and greater because of a 
partial saturation of the available number of activation centers in the 
material. Preflight calibration of the IMP detectors showed that this 
is best expressed as a function of the differential light output versus 
stopping power. In the energy interval of interest for this experiment, 
this functional dependence could be expressed as 

dL/dE = A - B/(dE/dx) (1), 

with B/A= 0.3. Integration of (1) in the case of the AE detector - 
produces the result 


AL 


A4E - B' 


( 2 ), 
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with the offset B' actually being a very slowly varying function of 
particle range and hence energy. Similarly a stopping particle of 
incident energy E - AE on the E - AE detector will give a light re- 
sponse of the form 


L = A(E - AE) - B' (3), 

in which again B* is a slowly varying function of energy. 

Because of the linear form of equations (2) and (3) with an offset 
it was possible then to express the overall response of the detectors 
A and B in the following form: 

E - AE = a + bx (4) 

and 

AE = f + gy (5). 

In equations (4) and (5) E - AE and AE represent the energy/nucleon 
deposited in the respective crystals, and x and y are the signal out- 
puts from the entire detector electronics train in terms of pulse 
height channel number. 

Since equations (4) and (5) each contain two undetermined con- 
stants, two calibration points along the response curve are required. 
The endpoint automatically provides the first of these calibration 
points. Ideally, the second point, or several points, is found by 
preflight calibration using a monoenergetic beam of particles. Such 
a calibration was carried out. The theoretical values for AE and 
E - AE for a particle with incident energy E were computed as de- 
scribed in Chapter II. These values of E - AE and AE corresponded 
to the channel numbers x and y found in calibration. This can be 
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located approximately by noting what channel the threshold of the E - AE 
signal occurs in. The corresponding AE channel is the y coordinate of 
the curve in the matrix having this x threshold value. More accurate 
determinations of several points along this distribution have also been 
made by utilizing a monoenergetic beam of protons from an accelerator. 

A combination of these two methods was used to calibrate both the scin- 
tillator and the solid state telescopes with agreement between the two 
methods. The coefficients for equations (4) and (5) for the various 
experiments are given in Table II. 

Using the Barkas and Berger (1964) range energy loss expressions, 
detailed tables were computed for the precise theoretical energy de- 
posited in each detector for a given incident energy for each one of 
the telescopes. Using these tables and the appropriate forms of equa- 
tions (4) and (5), it was then possible to ascertain the incident 
energies corresponding to a position on the matrix curves. 

Because of possible gain shifts in the satellite electronics, 
there could be a change in one, or both, of the pulse height analyzer 
responses. On IMP-III this was a very noticeable effect. In order to 
correct for such gain changes, the endpoint of every one orbit' matrix 
over the lifetime of the satellite was- examined and the E - AE and AE 
channel numbers of the endpoint determined in each case. Picking one 
such endpoint as the standard, all other endpoints, which always cor- 
respond to the same energy, could be compared to this standard and AE 
and E - AEgain factors could be computed. These gain factors could 
be applied to the data from the respective orbits, permitting combin- 
ing of the data for longer periods. Indeed this is the process that 
was used in combining the almost four months of data presented in Fig- 


ure 15 . 
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C. Matrix Histogram Analyses 

Consider again the region of the proton line in Figure 15. As 
noted, there are two outer curves bounding a region containing the 
proton line. This region was originally determined using only one 
orbit of data following a large solar flare so that the copious quan- 
tities of energetic protons so produced would easily delineate the line. 
It is readily apparent that there are a number of counts outside of this 
proton region. The origin of the non-proton events can be mostly attrib- 
uted to nuclear interactions in the telescope, and to a lesser degree in 
the spacecraft. These counts, though obeying the ABC criterion, are not 
primary protons and constitute noise that must be separated from the 
true primary events. An example of such a. reaction causing the back- 
ground counts could be a secondary neutron formed in the spacecraft by 
a cosmic ray primary. Such a neutron could penetrate the C detector, 
producing no signal in the process, and then decay in the B detector, 
where it would deposit a certain amount of energy and continue on 
through the A detector leaving there a small amount of energy. Such an 
event would certainly satisfy the ABC criterion, but would not be caused 
by a primary cosmic ray, which is what one is looking for. 

If one is careful to exclude regions containing protons, alphas, 
or electrons in the matrix, it is found that parallel to either axis in 
the matrix the "background" counts, as described above, obey a power 
law in channel number to a very good approximation. This background 
flux will contaminate the region of interest, and hence mu-st be removed. 
The power law behavior of the background provides a convenient method 
for removing this contribution in the neighborhood of the proton line, 
for instance, so that one can obtain a true measure of the primary proton 


counts . 
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In Figure 15 the region of the proton line has been divided up 
into six boxes which in Increasing order represent energy lines lying 
between the incident energy values of 18.7, 29.2, 39.7, 50.2, 60.6, 71.1, 
and 81.6 MeV. By taking narrow strips parallel to the proton line and 
lying inside of the upper and lower energy boundaries as shown, the 
number of counts in each strip can be determined. Extending the number 
of these strips in either direction beyond the proton region so that 
some of the background region will be included, it is possible then 
to plot histograms of the counts so determined. Such histograms cor- 
responding to the energy bins shown in Figure 15 are presented in Figure 
16 . The dashed lines show the mean values of the power law in channel 
number and were obtained by plotting the strip counts on log-log paper 
for several trial midpoints in channel number until the best straight 
line fit was obtained for the background regions. The shaded areas re- 
present the true primary proton counts in the energy interval specified. 
The region below the dashed line and between the vertical lines delimit- 
ing the proton region is the background component. 

As discussed in the previous chapter, the LED experiment was 
surrounded by a plastic scintillator cup constituting the C anticoinci- 
dence detector. Histogram analyses of the LED data analogous to that 
discussed above shows almost no background contamination. This justifies 
the assumption that most of the background secondaries were produced in 
the detectors, or the walls of the telescope housing, in the case of the 
scintillators by high energy primaries with kinetic energies > 1 GeV. 
Further, if histograms are produced for solar flare protons after a 
large flare, when very high fluxes of protons register as events, in 
the case of the scintillator telescope the relative background is 
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vanishingly small. This decrease in the number of secondary protons with 
respect to the low energy primaries is because the flux of secondaries 
has remained constant. This reinforces the view that secondaries are 
produced by nuclear interactions of high energy primaries which are not 
present in the solar flux and are not too highly modulated so that their 
flux through the spacecraft remains relatively constant. 

The differential proton flux was computed from the total proton 
counts in a specified energy box using the relation 


dJ i 


Mj N a 

N«j* t AE^ 


( 6 ) , 


where is the number of protons in box i, the total number of 
counts in the matrix, N a is the accumulated ABC counts during the 
accumulation time t, and and AE^ are the geometry factor and dif- 
ferential energy respectively for box i. The energy 'E^ taken for dJ^ 
is just the mean of the upper and lower energy bounds- of the ith box. 

Using this method of histogram analysis for both proton and alpha 
lines, it is possible to obtain very good counting statistics. In 
practice, preliminary examination of the data was made to determine 
the quietest time periods so that contamination by solar flare protons 
could be avoided. Excluding solar disturbed periods, data could be 
combined to cover several months such as the example shown in Figures 
15 and 16. For the large numbers of counts involved, Poisson statis- 
tics are applicable (See the discussion in Appendix C) . The relative 
error is a combination of the three relative Poisson standard deviation 
for the %, N a , and N T , i.e.. 


'_L JL -i 

% + N a + N t 


v 2 


cr 


(7). 
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Examination of the background shows that it is reasonably constant 
over the lifetime of the satellite. This fact provides a means of 
making a background correction to the short time period fluxes discussed 
in the next section. Such a background correction can be put in terms of 
a flux in the same units as the particle flux to be measured. This is 
accomplished simply by using in place of Nj_ in equation (6) the number 
of background counts below the respective dashed lines as shown in the 
histograms in Figure 16. This correction is to be subtracted from flux 
measurements in the same energy range for which the background counts 
have not been separated by means of the histogram method illustrated 
above . 

For the MED telescope the minimum flux measurable occurs when the 
flux and its computed standard deviation are approximately equal. 
Examination of 24 hour average proton flux values showed that this 

c o 

value occurred at about 10 _:, /s cm sr MeV. This is also about the 
average value of background flux that had to be subtracted from the 
MED fluxes to correct for nuclear interaction events in the spacecraft. 
Even though there was no background correction necessary for the LED 
experiment, the box method gave about the same minimum for a measurable 
flux. 

D. Computation of Particle Flux Time Histories 
With Fine Time Resolution 

The methods outlined in the last section for computing fluxes have 
the limitation that they are good only for time periods of the order of 
days to months. Since many of the cosmic ray phenomena of great interest 
occur in times of a few minutes to times of hours, some other method of 
computation must be used. In principle the same prescription as given 
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above can be followed, but since the size of the matrices involved is 
quite large, it is not feasible to construct these matrices every few 
minutes and then examine the histograms. Instead a computational 
procedure amenable to coding for a high speed digital computer was 
developed. 

The shortest time period records for the experiments were read 
from final data tapes and effectively the pulse height analyzed counts 
that fell into the various energy bin as depicted in Figure 15 were 
accumulated, as were the ABC counts, the total number of pulse height 
analyzed counts, and the ABC accumulation time. Two modes of compu- 
tation were available; the first being a fixed time period process in 
which counts were accumulated for a predetermined time and the fluxes 
and errors then computed on the basis of this accumulated data. The 
second mode computed the fluxes and standard deviations of the fluxes 
after the addition of each record of data, and continued to add more data 
until a predetermined precision criterion had been satisfied. These 
two methods were designated the fixed time and variable time techniques. 

To obviate constructing matrices of the ABC pulse height analyzed 
events as was done in the previously discussed histogram method, each 
event was examined to see where it would fall in such a matrix if it 
were to be constructed. This was accomplished by producing a mask of the 
region surrounding the particle line in question -for a standard line 
including the energy boundaries of the boxes as depicted in Figure 15. 

The mask was simply the x values of the box boundaries for" each value 
of y. For instance in Figure 15 if an event had the AE channel number 
of 40 and E - AE channel numbers of more than 31 and less than 38 then 
that event would be attributed to box 2. This mask data was tabulated 
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on a row by row basis for a standard line and used in the computer as 
a "go-no go" test in the manner described. 

One feature of this box counting of particle events was that it 
permitted the correction for gain changes in a fairly simple, manner . 

The gain factors as described in the last section were simply applied 
to the masks during the applicable periods. In this way there never 
was any need to consider fractional counts. 

Besides accumulating the individual box counts as described, the 
total number of ABC pulse height analyzed events was accumulated for a 
given energy interval being computed. Also the total ABC counts N g 
and corresponding counting time t were accumulated. It is fairly obvious 
that in principle the box counts obtained in this way and the total ana- 
lyzed ABC counts are equivalent to the total histogram particle counts 
and total matrix counts Nf, respectively, as descussed in the previous 
section. The only difference is that in the present case the background 
counts were included as well as the desired particles. To compute the 
particle flux using the method under discussion, equation (6) was used 
with replaced by the total box counts as described and Nrp by the 
total accumulated ABC pulse height analyzed counts. N and t were the 
same as before but accumulated over a very much shorter time period. 

The background rate as described in the last section was a constant 
throughout the lifetime of a particular satellite and hence was subtracted 
from the raw computed flux to give the absolute particle flux. For com- 
pleteness, then, the formula used to compute the flux by this method 
may be written as 

Nj N a 

N T t Gi AEj_ " r i 


dj 


( 8 ), 
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where rj[ represents the constant background component, and the other 
quantities are completely analogous to those used in equation (6) except 
for which contains background as well as primary counts. Also it 
should be noted that N^, although still the total number of ABC pulse 
height analyzed counts, is at most a small number depending, of course, 
on how long a period the averaging process extended for a single flux 
calculation. 

The use of Poisson statistics for this method of flux computation 
is not quite so clear cut. The ABC event counts were large and reason- 
ably unrestricted, hence the standard deviation in N a in (8) was given 
by the Poisson formula (See Appendix C, equation (35)) : 



and the relative standard deviation by 

S a = a a /N a (10). 

In the case of N^, however, the use of Poisson statistics is not justified 
because of the small number of bins N^, which could vary anywhere from 
4 to 8 per commutation cycle. As discussed in Appendix C, the appro- 
priate formula to use is derived directly from the binomial distribution. 
Using equation (19) from Appendix C, the standard deviation for is 
given by 


= [N £ (1 - %/Nt)]^ 

and the fractional standard deviation by 


Si = CTjVNi 


(ID, 


( 12 ). 
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The additional relation, equation (20) in Appendix C, was utilized where 
required. Treating the two errors as given by equations (10) and (12) 
as independent errors , the two may- be combined to give 

S = (S a 3 + Si 2 ) % (13), 

so that the standard deviation in dj may be written in terms of (13) as 

ct = S • dJ (14). 

The ’’box" method under discussion worked extremely well for obtain- 
ing proton fluxes since protons are the predominant component of cosmic 
rays-. For alpha particle fluxes which were also determined on a short 
time basis, one' drawback was a relative saturation of the electronics 
by high proton flux rates, thus making more difficult a true measure- 
ment of the alpha particle fluxes. One design feature in IMP- IV noted 
in Chapter II was a high threshold mode of operation which in principle 
would have suppressed a number of proton counts, thus enabling the 
counting of a higher number of alphas. As it turned out, the high 
threshold ABC sums were not transmitted from the satellite. Although 
an artificial ABC sum could be computed from the ratios of total pulse 
height counts in both modes and the low threshold ABC sum, the alpha 
pulse heights were still too scarce to compute reasonable alpha fluxes. 

The very scarcity of alpha particles noted made possible the use of 
a computational technique that utilized the number of alpha pulse heights 
in a given energy interval and the time interval that the pulse height 
bin was open to record an event. The derivation of this technique is 
given in Appendix D. Using equation (13) in Appendix D, the flux of a 
particle based on pulse heights alone Is given by 



44 


dJ = ,t®E ln CV ( % “ N i> ] (15), 

where t is the time the pulse height bin will accept a count, G is the 
geometry factor, and AE the energy interval. is again the number of 
pulse analyzed counts. Equation (15) is valid in the limit that « Nrp 
The standard deviation for dJ can be computed using equation (13) in 
Appendix D, so that 

CT = tGAE [Nj/N-p (Nrp - N^)]^ (16). 

One other statistical item was considered in the short period 
analysis. This was the problem of upper limits when there were no 
counts during the averaging interval. The procedure adopted when 
there was no count during an averaging period was to add one artifi- 
cial count and compute the appropriate flux on the basis of this one 
count. This flux was then considered to be the upper limit for that 
particular interval of elapsed time. 

In Table III is a summary of the analysis parameters -used for 
this study. These values were used for both the box method and the 
histogram method of analysis. 



IV. THE OBSERVATIONS 


A. Time Histories of Low Energy Proton Fluxes 
From June 1965 to August 1968 

The period concerned with in this study begins with the launch of 
IMP-III spacecraft near the minimum of solar activity in 1965. There 
is a certain amount of arbitrariness in speaking of a measure of solar 
activity. From a solar astronomer’s viexjpoint the minimum of solar 
activity for this current cycle occurred in October 1964 as indicated 
by the minimum in sunspot number shown in Figure 2. On the other hand, 
the cosmic ray flux at earth as indicated by neutron monitor counting 
rates reached a first maximum in April 1965 and a secondary maximum 
approximately six months later in October 1965. The current study is 
concerned with the behavior of the low energy component of cosmic rays 
during the period of decreasing galactic cosmic ray flux. 

Figure 17 shows a history of the 24 hour proton differential fluxes 
at energies of 29 and 60 MeV for the' entire period of this study. Using 
the results from the solid state telescope on IMP-IV a third flux value 
computed at 12 MeV is shown from late May 1967 up through the middle of 
August 1968 where this study ends. Because of the loss of communications 
on IMP-III there is missing data from 1 May to 24 May 1967. There is an 
additional block of data missing from thi-s analysis because of a bad 
data tape for the period from 10 April to 1 May 1967. This former is 
unfortunate for several reasons including the loss of data for the large 
proton event commencing on 24 May. 
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The gross features of the solar activity and production of ener- 
getic solar protons are summarized in Table IV. In this table a semi- 
annual count of the number of substantial increases in counting rate 
above the quiescent rate of protons is tabulated on an integral basis. 

The counts in each column are for fluxes greater than or equal to 
the column flux heading and therefore contain also the number in the 
adjacent column to the right. No account has been taken in this table 
as to whether the increase was flare related or a return to central meridian 
of an active region; however, the fluxes refer to the peak flux 
reached during any of the increases. The 29 MeV flux is derived from 
the lower two boxes of the MED experiment spanning the energy interval 
18.7 to 39.7 MeV hence, should certainly indicate particles of solar 
origin. 

Two aspects of the period of study are worth consideration. First, 
a maximum in activity as reflected by the high frequency of large events 
was reached in the first half of 1967. This was followed immediately 
by a minimum in number of events in the next six months and then a slow 
return to a higher number of events by. the end of 1968. This large 
maximum followed by the minimum in 1967 is probably related to a con- 
current decrease in quiet time fluxes to be discussed below. The second 
aspect of the frequencies shown in the table is the monotonically de- 
creasing behavior with respect to peak flux of the total numbers of 
counts over the entire period under discussion. This latter bears a 
resemblance to the frequency of optical solar flares as a function of 
flare area (Smith 1963). 

Returning again to Figure 17, it is seen that such events as are 
enumerated in Table IV represent enhancements of at least a factor of 
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10 above the quiet time rates. It is reasonable to assume then that 
because of the well known observed correlations of these short term 
proton flux increases with solar activity, either in the form of flares 
or returning solar active regions, that the origin of protons during 
such events is almost entirely solar. Conversely, it has been generally 
assumed that during very quiet solar periods when particle fluxes are 
at their minimum observed values, the composition of cosmic rays is 
predominantly galactic in origin. Such solar quiet times would be 
exemplified by the periods June through September 1965, October through 
December 1966, and most of September and October in 1967. It would be 
during these, or similar, quiet periods that one would want to study the 
long term variations of cosmic rays and in particular the solar modu- 
lations of the galactic cosmic rays. 

It should be noted here that all observations of solar flare 
associated proton spectra show a very steep monotonically decreasing 
behavior for proton fluxes as a function of kinetic energy. With this 
fact in mind and referring to Figure 1, it is readily apparent that 
the higher energy fluxes in this study should be better able to give 
good statistics on galactic cosmic rays during quiet times. On this 
basis the flux with the next to highest energy for these experiments 
and the one just below it are presented in Figure 18 on a monthly 
basis. The highest energy value was not used here because of some 
lack of confidence in its overall statistics. Plotted on the same graph 
are the monthly averaged hourly Deep River neutron counts. The proton 
flux averages are based on the quietest flux periods available each 
month. There are several months where no points are presented simply 
because of the greatly disturbed nature of that period. 



48 


In examining the cosmic ray quiet time behavior in Figure 18 
several marked features stand out. First of all there appears to be 
a lag of one to two months in corresponding features in the IMP proton 
fluxes with respect to the neutron rates. For instance, a big minimum 
in neutron counts occurs in August 1965, and the corresponding minimum 
in proton flux does not occur until October. The second maximum in 
neutron rates in November 1965 is followed by the corresponding proton 
maximum in January of 1966. In September of 1966 there is a massive 
Forbush decrease, which appears to be related to the abrupt drop in 
proton fluxes in October and November. In this respect it is to be 
noted that the proton fluxes do not recover to anywhere near their 
previous level. Referring back to Figure 17d, it is probable that 
this latter behavior is associated with the large solar flare occurring 
on 2 September 1966, with its resultant proton event and greatly in- 
creased proton fluxes that lasted the remainder of the month. This 
lack of recovery is possibly directly the result of a large change in 
the interplanetary medium following the September solar disturbance. 

It is not as clear cut as the lag just discussed, but there appears to 
be some small displacement in time between the two IMP proton fluxes. 

The observed lags seen here would tend to indicate a hysteresis effect 
in the proton fluxes with respect to energy, similar to that suggested by 
Balasubrahmanyan (1965), and Balasubrahmanyan, Hagge, and McDonald (1967). 

In Chapter V will be presented a model of the low energy cosmic rays 
based on the 4 day averages of differential fluxes spanning the energy 
interval from 4.2 to 81.7 MeV. This model and its analysis will be 
restricted to the IMP-IV data only which extends down to the lower ener- 
gies. A plot of these 96 hour averages for two of the lower energies 



49 


(one each from the LED and MED experiments) is shown in Figure 19. It 
is readily apparent that the fluxes for the two energies track each other 
quite well. Further it will be noted that the lowest- energy value varies 
considerably between widely separated extremes, as compared to the higher 
energy flux, on a percentage basis. This latter behavior may be inter- 
preted iri one of two ways. Either the lower energy particles are more 
highly modulated then at higher energies, or the lower energy fluxes 
represent a larger component of solar particles. Indeed this is the 
very question that the present study will attempt to answer. 

In view of the last point noted above, it is of interest to return 
to Figure 18 and consider the relative changes for the two proton fluxes 
represented there. It is easily seen that the short term month to month 
changes are larger for the higher energy flux. Further the large decrease 
in September 1967 produced a reduction in the higher energy flux of about 
2.5 times whereas the lower energy flux was reduced by a factor of less 
than 2.0. At face value this behavior seems to be just the opposite of 
that encountered above for the 96 hour proton fluxes at lower energies, 
in that the relative changes in fluxes is larger for the higher, rather 
than the lower energy. This is not necessarily the case, however. The 96 
hour averages must represent all of the data at the given energies with 
no selection of quiet periods, and also is at the low end of the energy 
interval where the spectrum has been observed to have a negative slope 
as shown in Figure 1. On the contrary, the monthly flux averages shown 
in Figure 18 were taken for the quietest periods available and lay well 
above the minimum appearing in Figure 1. 

The difference in the monthly averages of proton' fluxes at the two 
higher energies can also be satisfied by alternate explanations. If 
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the fluxes so presented are predominantly of galactic origin then the 
fact that higher energy flux changes proportionately more than the' lower 
must be interpreted as a greater effect of solar modulation on the higher 
energies in the limited range considered. If on the other hand there 
were a monotonically decreasing steady solar component inherent in these 
fluxes then the modulation of the lower energy particles could be the 
same, or greater than, the higher energy particles and this effect would 
be effectively masked because of the larger solar component at the lower 
energy. 

B. Proton and Alpha Particle Spectra 
And Cosmic Ray Modulation 

The time dependent behavior of cosmic rays is dominated by an 11 
year modulation that has been rather well established experimentally (see 
discussion in Chapter I) . By using the low energy MED fluxes it was 
possible to screen the data of both satellites for the quietest times 
available. Using this technique six periods as tested in Figure 20 ‘were 
abstracted from the data. Using the 'histogram technique described in 
the last chapter the proton and alpha particle fluxes were determined 
for these periods. The results are plotted in Figures 20a and 21a. 

The proton spectra presented in Figure 20a first of all show that of 
the data presented the June 1965 to October 1965 period and the November 
1965 to March 1968 period have about the same minimum amount of modu- 
lation. These two periods represent the maximum fluxes at the 1965 solar 
minimum and are in extremely good agreement with the same energy interval 
of the spectrum shown in Figure 1, which were derived by the Chicago 
group (Fan et al. 1966a) and the Goddard group (Balasubrahmanyan et al. 
1966a, b). The September to October 1967 spectrum is quite regular 
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except for the highest energy point, although two standard deviations 

X 

would allow this low flux. The most recent spectrum, the one. for April 
to May 1968, shows large modulation at the high energies but a rather 
peculiar double valued behavior at the low energy end. 

In Figure 21a again the June to October 1965 helium spectrum com- 
pared extremely well with the corresponding energy interval in Figure 
The fractional modulation between the protons and alphas comparing 
Figures 20b and 21b shows that on the energy per nucleon basis the pro- 
ton fluxes are reduced by a factor of 2.7 while the alphas are reduced 
by a factor of 3.5 during the same time. It should be noted that the 
November to December 1966 alpha spectrum shows, a very steep power law 
behavior as compared to the other alpha spectra. 

Using the formulation described in Chapter I in equations (3) through 
(7), the logarithms of the ratios of the minimum fluxes in 1965 to the 
fluxes for the other periods were plotted in Figures 22 and 23, for 
protons and alpha particles, respectively. As indicated by Webber (1968) 
the low energy cosmic ray rates should fall parallel to 1/3 functions. 

As seen in Figure 22 the proton logarithmic flux ratios have positive 
slopes unlike the negative slope expected. The alpha data in Figure 23 
shows even more erratic behavior except for the period of October to 
December 1966, and in this case the ratios are parallel to a 1/PP function 
rather than the 1/3 expected. 

In light of the evidence presented in the last section for some 
amount of hysteresis between fluxes of different energies, one might 
expect the erratic behavior observed. If such hysteresis does exist, 
then in order for the spectra and these modulations to return to the 
forms observed before solar minimum as described by 1|(t)/3 behavior. 
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then it would be necessary for a change in the formulation of the 
modulation. It is at least apparent that the accepted modulation 
theories do not apply during the period under discussion. 

As a tentative explanation to this problem consider equation (6) 
in Chapter I, 


dj p (P,t) = . AT] .. 

dj e (P,t*) f(P,P) 


(1), 


with ATI = Tl(t') - TKt). The exponent in equation (5) in Chapter 1 is 
just 


T K t ) = ^ w , 

f(P,P) J < dr (2)- 

1 A.U. 


Jokipii (1966) and Roelof (1967) have independently shown that the dif- 
fusion coefficient k in the integrand of (2) is related to the power 
spectrum of spatial irregularities in the interplanetary magnetic field. 
In other words, they have found a means of relating a property of the 
medium to the propagation of the particles. 

The form of the diffusion coefficient derived from theory may be 
expressed as 

- Pci 8 

3ttM(K q ) (3) > 


where pc is the .particle velocity, and K Q designates the wave number 
corresponding to a particle with gyro-radius r Q : 


K 


o 


1 




B 

P 


(4). 
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The power spectrum M(K) of the spatial variations is related to the 
time variations "seen." by a satellite by 

M(K) = £-m(f) (5), 

2tt 

where the frequency and wave number satisfy 

f « ~ K (6). 

2tt 

Measurements of the magnetic field in space (Webber, 1968), have shown 
that 

M(K) « 1/K, or 1/K? (7). 

With this form of power spectrum, the diffusion coefficient in equation 
(3) has the form 


< « PP, or .p (8), 

and hence in (1): 

f(P,P) = PP, or p (9). 

The contrary results exhibited by the proton spectra since solar 
minimum in Figure 22 show that the present theory must be modified so 
that 


f(P,p) « 1/P n 
or equivalently, from (7), 

1/ K nf2 


( 10 ), 


AM « 


( 11 ) 
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where the AM represents the change in the power spectrum since solar 
minimum until the measurement of the relative modulation, n must be 
1, or greater, then for the observed flux ratios to exhibit the form 
predicted by theory. 

C. Comparison of Quietest Time Energy Spectra 
From IMP- IV with Solar Minimum Spectrum 

Referring ahead to Figure 39, it is readily seen that during several 
time intervals in September and October 1967 the differential flux of 
protons at 5.2 MeV was at a minimum. Examination of these two months of 
proton and alpha particle fluxes in detail showed an interval from October 11 
through October 26 when the alpha particle counting rate almost vanished 
in the LED telescope. This period then was taken to represent the time 
of least solar activity on the basis of the low energy alpha particle 
fluxes . 

A combination of methods was used to obtain the fluxes during this 
period. Where counting rates almost vanished, but were still non-zero, 
the fluxes were computed from the few counts that fell into an energy bin 
and the total accumulation on-time for the period. In the cases where no 
counts occurred, one count was assumed. This one count was used to compute 
an upper limit to the flux. For higher counting rates, the box fluxes were 
averaged for the period in question. 

The proton energy spectrum obtained in the manner discussed is dis- 
played in Figure 20 (b), and the alpha particle spectrum is shown in 
Figure 21 (b) . Superimposed on each of these plots is a solid curve rep- 
resenting the corresponding 1965 quiet time composite spectrum (Gloeckler- 
and Jokipii, 1967). These composite curves are identical to the spectra 
indicated in Figure 1. As one would expect, the 1967 spectra occurring 
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during a time of- increasing solar activity show appreciable modulation. 
Also they show a low energy "tail" which, in light of the material to 
follow in Chapter V, is most probably of solar origin. If only the 
increasing parts of these 1967 spectra are compared with the 1965 spectra, 
one finds that the modulation of the protons and alphas is by a factor 
of 1.86 and 2.23, respectively, at 50 MeV/nucleon. 



V. A TWO COMPONENT MODEL FOR LOW ENERGY 


COSMIC RAYS AND ITS IMPLICATIONS 
A. -The Model 

Qualitative examination of the low energy proton fluxes during 
very quiet times in 1967, for instance in September in Figure 17 (g), 
show that there seems to be a very noticeable relative minimum at 
around 30 MeV. If a slightly more active period is examined, say in 
August 1967 in the same figure, it is apparent that the local minimum 
disappears as the lowest energy fluxes begin to dominate. Further, it 
is noticed that the lowest energy flux increases by several orders of 
magnitude while the adjacent flux increases by a much smaller amount 
and the highest energy flux by even less. This behavior immediately 
suggests that what one may be observing in the energy range considered 
is a monotonically decreasing spectrum of particles that undergoes large 
variations added to a relatively stabile monotonically increasing spec- 
trum of particles. The admixture of these components would seem to be 
such that during very quiet times the relative minimum noted falls in 
the vicinity of 30 MeV. It is probably not fortuitous that the proton 
spectrum in Figure 1 exhibits a minimum at about this same energy. 

It is very difficult not to immediately identify the decreasing low 
energy component as originating in the sun and the increasing component 
as being of galactic origin. Solar proton fluxes that are known to be 
associated with flares all have very steep monotonically decreasing 
spectra. Further, it has been generally assumed that the higher energy 
particles were of galactic origin as in Figure 1. It does not seem 
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unreasonable to assume then that at least part of the cosmic rays below 
the relative minimum in Figure 1 are of solar origin. Fan et al. (1968, 
1969) have suggested that the lowest energy proton fluxes observed by 
them below 10 MeV might contain a solar component. 

With the above observations in mind, it is then possible 1 to con- 
struct a two component model of the low energy cosmic rays which , as 
will be seen later, fits the observations remarkably well. Solar flare 
proton spectra have been shown to fit an exponential rigidity law for 
a wide range of rigidities (Fichtel and McDonald, 1967), but over a 
limited energy range can equally well be fitted to a power law in 
energy, e.g. 

J(E) = J 0 E _S (1). 

It would seem reasonable then to assume the same form for the quiet 
time solar fluxes that would apply to the model being discussed. 

Further, a power law in energy with a positive exponent could equally 
well apply to the increasing, or galactic, component in the limited 
energy interval being considered, even though above this interval it 
will no longer hold true. Very simply, then, such a model may be 
expressed as 

J(E) = F g E _ S + F q E G (2), 

where Fg and Fg are the values of the fluxes at unit energy, for the 
solar and galactic components, respectively, and S and G are the 
magnitudes of the power law indices in the two cases. Differentiating 
(2) with respect to E and setting the result equal to zero, one ob- 
tains for the energy at minimum flux 
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E •„ - (SV<5F G ) 1/(< *- S) 


(3). 


Asympo tic ally on either side of the minimum, the flux becomes dominated 
by one component, or the other, since 


and 


J 



for E < E . 

mm 



for E > E . 

mm 


(4), 


(5). 


The four independent parameters describing this model, Fg, Fg, S, and 
G, all may be functions of time. It is to be expected that Fg and S will 
vary more rapidly since they are related to transient solar phenomena, 
and the galactic component will show on the average only a comparatively 
slow variation, provided the averaging time is long enough to remove 
transient modulation effects. If indeed the galactic component is 
relatively constant for short periods, then the energy at which the 
minimum 'is observed should be almost entirely dominated by the varia- 
tions of the solar component. Eliminating Fg between equations (2) and 
(3) one obtains the result that 


^min 


(1 + G/S) F E E 


G 

min 


( 6 ). 


If the galactic component is represented by Jg then equation (6) may be 
expressed as 


J G * 1®< S+ G)]J ml „ (7). 

Hence, if indeed the observed spectrum is due to a rapidly varying solar 
component plus a relatively constant galactic component, then equation (7) 
is to be taken as a prescription for computing the galactic component 
from a series of measurements of J min , Emin> S ’ and G ' 
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It is now possible to summarize the assumptions and predictions of 
this model. With only the basic assumption that equation (2) represents 
the cosmic ray fluxes in some limited energy range containing the rela- 
tive minimum it follows that 

1. The lowest-energy fluxes and highest- energy fluxes should be 
relatively independent of each other, 

2. On either side of the relative minimum, the fluxes should 
asymptotically approach the simple power- law form as given by 
equations (4) and (5). 

If in addition it is also assumed that the solar component will exhibit 

> 

greater magnitude and more rapidly varying time dependence than the 
galactic component, then the following are also predicted: 

3. The location of the relative minimum will be closely related 
to the lowest energy fluxes, and almost independent of the fluxes 
at high energies, as long as that minimum falls between the fluxes 
at the two extremes. It should be noted that when the solar com- 
ponent becomes large enough, it should even dominate the higher 
ene rgy f luxes . 

4. The envelope of the various minimum fluxes taken over a period 
of time should trace out a power *law function related to the galac- 
tic component by the relation given in equation (7). 

With the above picture of the two component model in mind, it was 
possible to compare this model to observed proton and alpha particle spec 
tra. 

B. Least Squares Fitting of Model to Data 
In order to compare the model presented in the last section to the 
observed data, it was necessary to systematically fit the mathematical 
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formulation given to this data in order to determine the four parameters 
Fg , Fq, S, and G. Several independent methods of curve fitting were 
tried out before one method was adapted as giving the most self con- 
sistent results. The method used was conveniently available in the 
form of a Fortran language subroutine that made it possible to perform 
the calculations on a high speed digital computer. 

The algorithm used is completely general in that no special func- 
tions must be used to express the desired function. The method is fur- 
ther very general in that any function which is a good measure of 
goodness of fit between the computed and observed values, such as "chi- 
squared," may be used for computation. The process in essence steps 
each parameter of the desired function iteratively until the fit mea- 
sure reaches a relative minimum. After variation of one parameter 
achieves a minimum, then the next parameter is varied. This procedure 
is carried through for two cycles of variation using a process of para- 
bolic interpolation about the minima and increasing and decreasing the 
step size as necessary. 

For convenience it was decided to use as the measure of goodness 
of fit an expression for "chi-squared" given by 

N 

X 3 = 2 Cf(xi) - yiP/dyj 3 (8), 

i=l 

where fCx^ represents the theoretical function to be fitted, the y. are 
the dependent observed variables, the x^ are the independent observed 
variables, and the dy^ are the standard deviations or measured errors. 

N is the number of data points so that if P represents the number of 
free parameters in f (x^) , then the distribution represented by x 3 should 
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have a maximum at 


X 2 = N - P (9). 

By using x 3 then as the function to be minimized, it was possible to 
obtain at the same time a direct measure of the goodness of fit. Ex- 
plicitly 'the form of equation used was 

10 g 

X 2 = S [F s E i " + F g E i G - J i ] 3 /(7 i 2 (10). 

i=l L 

Ideally it would have been of great interest to treat as short time 
averages of the fluxes as possible in order to ascertain a high resolution 
picture of the time dependence of the model. In an initial attempt 24 
hour averages of the 10 differential fluxes for the protons and 7 values 
for the alpha particles were tried, but it was found that the overall 
statistical accuracy was too poor to give meaningful results. It turned 
out that 96 hour averages of the same data gave very good results and 
hence the- resulting analysis of the data with respect to the model was 
carried out on this basis. 

Since the computations were performed on a digital computer, it was 
possible to code into the same program the necessary output formats to 
make interpretations as simple as possible. The analysis was applied 
only to the IMP- IV data covering the period from 24 May 1967 to 18 
August 1968, so that the model could be examined at as high a resolution 
as possible in the energy interval from 4 to 80 MeV. Since there were 
451 days of continuous data included in the analysis, there were 113 
separate spectra each for the proton and alpha particle fluxes that were 
fitted. It is not particularly illuminating to merely display 113 graphs 
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of the fitted curves, or alternatively 113 sets of the four parameters 
Fg, F g , S, and G, although .these were produced directly and have a 
certain value. Rather it was expedient to perform directly on the 
fitted spectra the requisite comparisons and additional calculations 
to concisely compare the results with the predictions of the model. 

The procedure was as follows. First the least squares fitting as 
described was carried out providing the best values of the four para- 
meters in question for each of the 113 proton and alpha particle spectra. 
Next the minimum flux and the energy at minimum flux were determined 
numerically using the fitted values of the function. For each 96 hour 
set of proton and alpha spectra, a table was produced containing among 
other quantities the observed fluxes, energies, standard deviations, 
computed fluxes, final values of the function parameters, minimum flux, 
energy at minimum flux, and x 2 • At the same time a graph of the com- 
puted spectra for both the protons and alpha particles was made. It 
contained the data points with error bars superimposed on the fitted 
curves. In addition the ratios of helium to proton fluxes were computed 
for a number of values of energy/nucleon and a graph of these results 
were produced. 

After the above computations were performed, the results were 
treated collectively in order to examine the multitude of computed 
results. Graphs for both protons and helium nuclei were produced show- 
ing the following correlations: 

1. The flux at maximum energy versus the flux at minimum energy. 

2. The flux at minimum energy versus the energy at minimum flux. 

3. The flux at maximum energy versus the energy at minimum flux. 

The minimum flux versus the energy at minimum flux. 


4 . 
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5. The power law indices versus the flux at minimum energy. 

6. The power law indices versus the flux at maximum energy. 

With the results in the above described form, it was then possible to 
examine the results in a meaningful manner. 

C. Interpretation of Results 

Since it is not feasible to present all of the graphic results of the 
curve fitting herein, four representative graphs showing proton and alpha 
particle spectra are shown in Figure 25. Going from left to right and 
top to bottom, these show the actual data and computed spectra for four 
periods having progressively less amounts of the low energy components. 

As would be expected, these correspond to progressively quieter times of 
solar activity. It will be noted that indeed these log-log plots of 
observed data show straight line behavior on either side of the mini- 
mum. Further, it is to be noted that the high energy component is com- 
pletely masked by the high level of the solar component for the 13-17 June 
period. Also in the other three graphs the level and slope of the galactic 
component remains quiet constant. This latter is consistent with the 
assumptions made for the model. The alpha particle spectra are seen to 
have essentially the same behavior as those for the protons. 

Figure 26(a) displays four representative alpha-to-proton ratio plots 
for the same dates as Figure 25. All of the ratio plots show large time 
variations with no apparent regularities. A plot of the means and standard 
deviations for several discrete energies is shown in Figure 26(b). Two 
important features exhibited by this long term average are 

(1) A maximum ratio of about .2 at 50 MeV/nucleon, 

(2) A decrease in the size of the variations for increasing energy. 
Studies of flare associated events (Fichtel and McDonald, 1967) show that 
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the alpha-to-proton ratio varies between about 10"^ and lO"- 1 -. At galactic 
cosmic ray energies (Webber, 1967) the ratio is relatively constant at 10"^ 
after falling from about .2 at 60 MeV/nucleon. The decrease in size of the 
fluctuations with increasing energy may be interpreted as the lower energy 
ratios being primarily due to solar particles, while the higher energy ratio 
is more likely of galactic origin. The observed maximum in the ratio and 
decrease above 50 MeV/nucleon is consistent with observed ratios for 
galactic particles. 

As has already been suggested, the lowest energy flux should prob- 
ably consist mostly of solar particles and the highest energy flux of 
galactic protons. The 76.4 MeV proton flux is plotted versus the 5.2 MeV 
flux in Figure 27. Up to about a flux of 10"^ protons/s cm^ sr MeV at 
5.2 MeV, there appears to be no correlation between the two energies. 

Above this flux value the correlation becomes quite good. This behav- 
ior is consistent with the prediction that during quiet times the two 
extremes of the spectrum come from independent sources. 

With the assumptions made in the first section, the model predicts 
that for a relatively constant galactic component that the lowest energy 
flux should correlate very well with the energy at minimum flux while the 
flux at the highest energy would show poor correlation with this energy. 

This is very strikingly the effect that is brought out in Figures 28 and 
29. Figure 28 shows extremely good correlation between the 5.2 MeV flux 
and the energy at minimum flux; while the following figure shows almost 
no correlation between the 76.4 MeV flux and the energy at .minimum flux. 

In Figure 30 is shown a monotonically increasing dependence of the 
minimum flux as a function of the energy at minimum flux. As pointed 
out earlier, if the model is correct then such an increasing dependence 
should be related to the galactic component by equation (7). The slope 
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of the distribution in this plot is approximately 1 which is in extremely 

good agreement with the slopes of the two quiet time galactic components 

-4 2 

shown in Figure 25. If one takes a value of J s 10 at 10 MeV in 

min 

Figure 30, and nominal values of G and S of 1 and 3 respectively, then 
from equation (7) J^, — .75 J m ^ n — *75 x 10“^ which is quite close to the 
values found in the last two plots in Figure 25. 

In the next two plots. Figures 31 and 32, are shown the distribu- 
tion of the power law indices as a function of the low energy flux and 
high energy flux, respectively. With respect to low energy flux, there 
appears to be no systematic feature in either index except for an indica- 
tion of a very slight anticorrelation between the positive values of S 
and the low energy flux. This latter can be interpreted as a hardening 
of the spectrum during times of higher solar activity. The correlations 
of the indices with high energy flux are either non-existent, or very 
small . 

It would be expected from the behavior exhibited qualitatively by 
the alpha particle spectra. in Figure 25 that essentially the same behav- 
ior would hold for these particles as well. Indeed this is so for the 

most part with a few significant differences. Figure 33 shows no correla- 

_3 

tion between the high and low energy flux up to about a value of 10 
particles/s cm^ sr MeV/nucleon for the 5.2 MeV flux. Above this value 
there appears to be a small amount of anticorrelation which could be 
interpreted as the onset of modulation of the higher rigidity alphas 
because of the greater solar activity. Figures 34 and 35 exhibit the 
comparatively much greater activity of the solar alphas as compared to 
their galactic counterparts. The very flat spectra for the quiet time 
alphas shown in Figure 25 is certainly well supported by the flatness 
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of minimum flux distribution with respect to energy presented in Figure 26. 
Finally, there is no indication of any significant regularities in the 
distribution of the power law indices of the alpha particle spectra with 
respect to either the low or high energy fluxes. 

D. The Relation of the Two Component 
Model to Other Observations 

In the last section it was shown that the two component model is 
quite capable of describing the observations of proton and alpha fluxes 
in the vicinity of the relative minimum in the cosmic ray energy spectrum. 
Whether this model is justified on other than the phenomenological grounds 
of the data fitting the functional form presented depends primarily on 
some substantiating evidence that the sun is indeed producing particles 
with energies of at least up to about 30 MeV and on a quasi-continuous 
basis. From the very high correlations between large solar flares and 
the subsequent arrival of large, proton fluxes, there is little doubt 
that the sun produces such particles on occasion. In the absence of 
flares, however, the question arises of what, if any, features and prop- 
erties are present which would oermit a close association between ob- 
served enhancements in particle fluxes and the sun as their source. 

Examination of any time history of cosmic ray fluxes which is taken 
at low energies and with averaging times of from a few minutes to a few 
days will show that there are many flux increases of varying sizes and 
shapes. A few of these with very fast rise times, large peak fluxes, 
and exponential decays, usually can be immediately associated with 
solar flares with significant importances. The remaining events must 
be scrutinized in some detail, testing for periodicity' and correlation 
with solar features other than flares. This latter process usually 




67 


permits assignment of the bulk of the non-flare events to either 27 day 
recurrence type flux increases, or less often to isolated events show- 
ing no tendency to reappear. These latter two classes of events are 
certainly closjely associated with magnetically active regions in the 
solar atmospheres, which may be responsible for both the acceleration 
and long term storage of the particles. 

In Figure 39 are presented the 24 hour averages of proton fluxes 
at 5.2 MeV over the period of analysis for IMP-IV. Events associated 
with confirmed flares of 2N to 2B importance, or greater, are marked 
with an F. Most of the remaining events have been associated with 
particular large area calcium plage regions whose central meridian 
passage on the sun nearly coincides with the event. The plage region 
numbers are used to label each of these events in the figure.. The strik- 
ing feature of these events so identified is their close association with 
the return of the same region for several solar rotations whi ch are indi- 
cated by connection of the numbers of the returning regions with hori- 
zontal broken lines. On about 15 June and 9 July two isolated events 
which are not recurrent are indicated. 

A better idea of the character of these two events is shown in 
Figure 40. The one shown in June is seen to be very symmetrical with 
slow rise and decay times on either side of a rather constant plateau 
at the lower energy. This is certainly the behavior one might expect 
from a broad stream of particles of constant intensity rotating past 
the observer. In contrast to this event, a small flare event occurred 
on 3 June and shows the characteristic rapid rise time and slow decay. • 
The other isolated event occurring about 9 July is not nearly so -sym- 
metrical as the one in June, but again shows the ’slow rise and fall 
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with an intervening plateau. About 2 August a rather irregular prolonged 
flux increase occurs which is identified in Figure 39 as a 27 day recur- 
rence event. 

The IMP- IV data just presented shows considerable activity during 
the period considered. Although the two component model was applied only 
to this data, at a time when solar activity was quite high, it is of 
interest to see if at a quieter time the sun could also be a potential 
source. With this end in view, the integral flux of protons with energies 
greater than 20 MeV obtained from IMP-III is plotted on a six hourly basis 
from March 1966 to March 1967 in Figure 41. As in Figure 39 the flares 
and calcium plage region meridian transits are marked where these can be 
shown to be related to the observed flux enhancements. It is evident 
that even during the relatively quieter times .represented there are re- 
currence events that are periodically adding their output to the particle 
populations . 

Table IV showed that the number of events increased quite markedly 
the smaller the size of the event considered. From this one might con- 
jecture that there may be even more events below the level of detection. 
Such a conjecture is certainly not necessary in the case of the IMP-IV 
data t-o which the two component model was fitted. For the low energy 
flux data available there is an ample supply of solar protons, as shown 
in Figure 39, to keep the low energy component of energy spectra at an 
appreciable level. It will be necessary to make measurements during much 
quieter times, such as prevailed just after IMP-III was launched and to 
see how these fit the model- in order to determine if the low energy "tail” 
of the spectrum is greatly reduced. 
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There exists at least one case in the literature in which proton 
flux observations made at a different phase of the solar cycle seem to 
be compat able with the model up to somewhat higher energies. Meyer and 
Vogt (1963) present a series of differential- proton energy spectra 
obtained from balloon-borne detectors flown in July and August of 1961. 
They contend that these spectra which are very steep monotonically de- 
creasing functions of energy are primary cosmic rays. The relative 
minima exhibited in four of these spectra occur from about 180 MeV 
up to about 300 MeV, which is much higher than the minima observed 
since 1963 which have been in the vicinity of 30 MeV. Further, the 
lower energy decreasing part of their spectra which extends in energy 
from 100 MeV to the minima exceed by about an order magnitude the gen- 
eral flux levels observed from 1963 to the present time. 

The above spectra of Meyer and Vogt are much more easily explained 
on the basis of the two component model' presented here than by assuming 
that they are indeed galactic primaries. If the model is correct then 
they observed the high energy end of the solar flux component where it 
was of comparable value to the primary component resulting in the rela- 
tive minimum. This interpretation of their results does not disagree 
with the IMP- IV results, since for such solar component levels the mini- 
mum would have fallen above the IMP-IV energy interval. 

In closing it should be emphasized that the two component model 
presented is by no means confirmed by what has been discussed. Further, 
it has not conclusively been shown that the lowest energy fluxes are all 
solar protons. Rather with respect to the first of these items, it 
would probably suffice to say that the necessity condition of a proof 
has been met, but not the sufficiency criterion. Much further testing 
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of such a model will be required over broader energy intervals and for 
the entirety of the solar cycle. As to the second point, it is prob- 
ably clear that any increases above the quiet time flux levels are 
of solar origin. When this latter is accounted for, what remains 
is asserted to be of galactic origin. 



APPENDIX A 


GEOMETRY FACTOR FOR A TELESCOPE WITH CYLINDRICAL SYMMETRY 

Consider the schematic of a particle telescope as shown in Figure 6, 
of height h and base radius Ri and top radius Ra Assume that the 
particles have an isotropic flux of J particles per unit area and solid 
angle. Referring to the diagram, then the differential number falling 
on the increment of area dA and passing through the upper end area 
increment dA' will be given by 

dN = J cos 9 dA dr (1) 


where 0 is the ‘angle between the vertical h and s, the line connecting 
the two area elements, and where 


dA = rdrd0 
dA 1 cos 9 


dr = 


and 


dA' = r'dr ' d0 


From the diagram it is also seen that 


l 2 = r^ + r' 2 - 2rr' cos(0* - 0) 


s 2 = h 2 + l 2 


and 


cos 9 — h/s 


( 2 ) 

(3), 


(4). 


(5) , 

( 6 ) , 


( 7 ). 


Equation (1) can now be integrated using these relations with the 
result that 
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N = I G. (8) 

where the geometry factor G is 

Ri Eg 2 tt ?H-tt 

G = 2h^ J rdr J r'dr' J d0 J* [h^ + r^ + r'^ - 2rr’ cos (0 1 - 0)] - ^ d0' (9) 
0 0 0 0 

The four integrals in (9) are easily performed going from right to left 
with the final result 


rfh 3 

r, , Ri 3 + fe- 2 

4. 0 Rl8 + ^ J. (E ’ 3 - Ep 3 ) 3 ] 

y a ) 

2 

i 1 * ■ 

l h 3 h 4 J 

2 1 


( 10 ) 


If the linear dimensions are measured in centimeters and the angles in 
radians, the units of G are centimeters 3 steradians (cn? sr) . In the 
case of a right circular cylinder of radius R, i.e. Rj. = Rg = R, equation. 
(10) reduces to the form 


rfh 2 i 

f 2B? 

f 41?) / 

£ " 

G = 2 1 

1 + * - 

M 



( 11 ) 


For a long slender telescope, i.e. where R/h « 1, equarion (11) reduces 
to the approximation 


G 


if R 4 /h? 


( 12 ). 
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APPENDIX B 

COMPUTATION OF AVERAGE PARTICLE PATH LENGTH 
IN CYLINDRICAL GEOMETRY 

If the line s in Figure 6 represents a typical particle track for 
a particle making an angle 8 with the axis, then to find the average 
-track length <s> for an isotropic distribution of particles it is 
necessary to average all such tracks over an appropriate distribution. 
This is very easily done if it is noted that in Appendix A the compu- 
tation of the geometry factor G is accomplished by integrating over all 
angles within the limits such that the track s passes through the two 
ends of the cylinder. 

Using Appendix A we can write for G 

h 3 dA dA' n \ 



tirhere s =; s (r , r ' ,0,0' ) . It is then possible to write a normalized 
distribution function for the angles in the form 


w(r,r' ,0,0’) 


h 3 

G s 4 (r ,r ' ,0,0' ) 


( 2 ). 


The average of s over an isotropic distribution is given then by 


<s> = w(r,r* ,0,0' ) s(r,r',0,0') dA dA' 

h 3 pp dA dA' 


<s> 




(3) , 

(4) . 


Using the equations (2) through (7) from Appendix A this can be written 


in the explicit form 
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% Eg 2n 0+TT 

I rdr J* r ' dr ' I d ^ J ■ 

0 0 0 0 1 


[h 2 + r 3 + r 13 - 2rr'cos (0' - 0)] (5) 


In theory equation (5) is integrable. It was found in practice that it 
was more convenient to perform the integration numerically on a digital 
computer over the range of required values. This was done quite easily 
using a Monte Carlo integration algorithm. 

The computed values were expressed as an average value of the cosine 
of the angle between the average path length and the longitudinal axis of 
the telescope such that 

<cos 0> = h/<s> . (6). 
The corrections for both telescopes in this study were only a few percent: 


1 / <c os 9> — 1.03 


(7). 



APPENDIX C 


THE BINOMIAL AND POISSON DISTRIBUTIONS 


For most random processes, including the counting of cosmic ray 
partic les Poisson statistics are usually assumed. Since the statis- 
tical requirements of the several forms of data analysis described in 
the text are varied, it is pertinent to derive here the requisite for- 
mulae for computing the standard deviations of these experiments. 

One is concerned with the counting rate r which is by definition 
the number of particles n counted in the time interval t, i.e. 

r = n/t (1) 

If the interval t is divided into b bins such that the probability of 
any two counts occurring in one bin is negligible, then one can write 
the probability that a count will occur in a bin as 


P = n/b (2) 

where n is the average number of counts in the time t. Similarly one 
can write the probability that such a count will not occur in a bin as 


q = 1 - p = 1- n/b 


(3) 


The combined probability that n counts will occur in the time t, or b 
bins, in a particular order is simply 


n _ 

fa] 

n 

h - *1 


[bj 


l bj 


b- 

(4) 


Now there are b! different ways that such a particular order can be 
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obtained for distinguishable counts. There are n! different ways that 
the n counts may be arranged among themselves and (b - n) ! different 
ways the remaining bins may be ordered among themselves. It follows 
then that the number of ways n indistinguishable counts can occupy b 
bins is given by the expression 

b_I 

n! (b - n)! 



From equations (2) and (3) 


p + q “ 1 

so that 


b 

E w b (n) = (p + q) n = 1 

n=0 


(7), 


( 8 ), 


and hence w^ (n) is the probability distribution for counts occurring 
in b bins. By definition the mean value of n counts occurring in b bins 
is 

b 

n = E n w fc (n) (9). 

n-0 
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Using equation (6) one may write the expression 

(py + q) b = 2 1^] p n y n q b ~ n = 2 y n w b (n) 

n=0 ^ n=0 b 


Differentiating equation (10) with respect to y one obtains 


- (py + q) b = b(py + q) b 1 p = 2 nw, '(n) y n_1 

d y n=0 


Setting y equal to 1 in equation (11) the expected result is 


n = 2 nw (n) = bp 

n=0 b 


Differentiating equation (11) again with respect to y one obtains 


b(b - 1) (py + q) b ~ 2 p 2 = 2 n(n - 1) w fo (n) y 11 2 

n=0 


Setting y equal to 1 in equation (13) the result is by definition 


mean value of n(n - 1) , or 
b 

n(n - 1) =2 n(n - 1) w^ (n) = b(b - l)p 2 

n=0 


Expanding the left hand side of equation (14) one obtains 
n(n - 1) = n 3 - n = b(b - l)j? 

By definition the variance is 

= (n - n)^ = n 2 - ~x\ 

where cr is the standard deviation. Substitution of equation (12) 
(15) into' (16) produces the result 


( 10 ). 


( 11 ). 


( 12 ). 


(13). 

the 


(14). 


(15). 


( 16 ), 

and 
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2 

a 


bp(l - p) 


(17), 


so that the standard deviation for n counts .in b bins is simply 

CT n = [bp(l - p)]Vfe (18), 

or using the probability of a single count occurring in a bin as given 
by equation (2), (18) becomes 

CT n = [nd-n/b)] 1 ^ (19). 


Equation (19) is a completely general formula for the standard 
deviation of the random process of putting n items in b boxes. As 
will be shorn shortly, the Poisson distribution and its resultant sta- 
tistical quantities are special cases following from the binomial dis- 
tribution given in equation (5), with the condition that b and n are 
both very large. The binomial distribution, on the other hand, is 
applicable in the general instance and hence holds for small b and_ 
n. One inadequacy in the binomial standard deviation as given by 
equation (19) is that cr n vanishes for n = 0 and n = b. In practice 
the standard deviations may be given realistic values for these two 
cases by arbitrarily adding one count and one bin to n and b, respec- 
tively. Equation (19) is in these two cases 

cr n = (n+ 1) [i- (n+ l)/(b+ 1)] Vz (20). 


If one considers the case where b becomes very large then one can 
take the limit of equation (5) as b approaches infinity. In this limit 
the binomial coefficient becomes 



b(b - 1) (b - 2) . . . (b - n + 1) = b^ 


lim 

b -4- co 


lim 
b -»• 03 


( 21 ), 
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and the right most factor in -equation (5) is 


lim f 

1 - s 

| b ' n = lim 

[l - 

* 

(b-ri) U 

+ ‘(b-n) (b-n-1) 

n 3 _ 

,n b -n' 

I 

|j CO l 

bJ 

b-» 03 


1! b 

2! 

b 2 

* b b - n . 


= e" n (22). 

With the two limits obtained in equations (21) and (22), the limiting 
case for equation (5) is given by 

w(n) = lim (n) = n 11 e n /n! (23), 

b -M» 

which is the Poisson distribution function. 

In the limit of Sterling's approximation which is valid for n 2 10 

n! £s /2 tt n e (24) , 

one may express the logarithm of (23) as 

In w(n) - n In n - n - In /2 tt - (n + 1 / 4) In n - n (25). 

Differentiating equation (25) with respect to n and setting the result 
to zero one obtains 

lnn-lnn-(n+ 1 / 4)-i+l = 0 (26), 

n 

or 

n = ne^ n (27). 

It is easily seen that in the limit of large n, 

n — n (28). 

In this case of large n, if one expands In w(n) in a Taylor's series 
about n, the first derivative with respect to n vanishes because of (28) . 
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Keeping only the first non- vanishing terms above the constant term, one 
obtains 


In w(n) = 


In w(n) + 


[ cP In w(n) ') 

dip 


X 3 


(29) 


where x = n - n. Using Sterling's approximation for n, the first term 
in (29) becomes 

1 


w(n) = 


n n e -n 


n -n 
n e 


n 


i / 2rm n“ n e n /2 tt h 


(30). 


Differentiating equation (25) twice with respect to n, the result is 


d 3 In w(n) _ 1 _ 1 

dn 3 2n 3 n 


(31) 


Substituting equations (30) and (31) into (29) one obtains the result 

In w(n) = - In /2rr h - — j — -—lx? (32), 

2 l 2n 3 H J 


or keeping only the largest second order term this may be. expressed as 

(33). 


I -x 3 /2n 
w(n) = — ^=== e 


/2nn 


Comparing equation (33) with the usual Gaussian error distribution 


£(x) = 575S 6 


-x 3 /2 CT 3 


(34), 


it is easily seen that the Poisson standard deviation is just 

a = /n~ (35). 

Since long counting runs are usually encountered where Poisson 
statistics are valid, it is the usual practice to replace n by n in 
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(35). It should be noted by comparing equation (19) and (35) that for 
small b one cannot use Poisson statistics as is often erroneously done. 
As one would expect, however, equation (19) in the limit of large b 
reduces to equation (35) . 



APPENDIX D 


COMPUTATION OF COUNTING RATES 
BY THE DEAD TIME APPROXIMATION 

In the case a time bin is open for a fixed period of time to accept 

one, and only one, count during this period, it is possible to compute 

a counting rate if the probability of a count occurring is much less 
than unity. Let r be the counting rate that is sought, and let t be the 
period such that a bin may be open to receive such a count. Let P(l,t) 
be the probability that one count will occur in the interval after the 
time t. Then this may be expressed as 

p(i,t) = p ( 0 , t) • P(i,t,dt) (1), 

where P(0,t) is the probability that no count occurred before the time 
t, and P(l,t,dt) is the probability a count occurred in the differential 
time dt at t. If r is the counting rate then the second factor on the 
right hand side of (1) may be expressed as 

P(l,t,dt) = rdt (2). 

If the interval t is broken up into N intervals dt = t/N, then the prob- 

ability a count does not occur in t is given by 

P(Ojt) = J Pi (0,t i ,dt i ) (3), 

i=l 

where the probability that a count does not occur in dt^ and t^ is -just 

P £ (0,t.,d ti ) = 1 - rdt. (4). 
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Using equation (4) in (3), one may write 


P(0,t) = J (X - rdt . ) = lim (1 - rdt) N 

17 1 N -» co 

i= 1 


lim (1 - 
N-* 00 


= e 


•rt 


(5) 


Using the results of equations (2) and (5) in (1), one finds that 
the probability for one count occurring after a time t is just 

P(l,t) = re' rt dt (6) 

The total probability that one count will occur sometime during the 

interval t is by integrating (6) over this interval, 

t 

P(t) = J* re“ rt ’ dt' = 1 - e“ rt (7) 

0 

If there are b time bins of duration t and n counts such that n « b, 
it is also true that 

P(t) = n/b (8) 

Combining equations (7) and (8), one obtains 

e~ rt = 1 - n/b (9) 

or 

r = — In [b/b - n)] (10) 

for the counting rate. 

If equation (10) is differentiated with respect to n, the result 

1 

t(b - n) 


then 


dr = 
dn 


( 11 ), 
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is obtained. The differential may be replaced with the errors, or 
standard deviations in this case, of r and n, hence 


gn 

a r t(b - n) 


( 12 ). 


CT n can be evaluated using equation (19) in Appendix C, with n replacing 
n, so that 


_ I n x / z - 
t b(b - n) 


( 13 ). 
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TABLE I 

SATELLITE DATA 



• IMP- III 

IMP- IV 

Launch Vehicle: 

3-stage Delta Vehicle 

3-stage Delta Vehic 


DSV-3C 

DSV-3E 

Launch Site: 

Eastern Test Range 

Western Test Range 

Launch Time: 

0700:00 May 29, 1965 

1405:54 May 24, 196' 

Weight: 

59.4 Kg 

73.8 Kg 

Initial Perigee: 

208 Km 

248 Km 

Apogee: 

260,800 Km 

2 1 1 , 1 16 Km 

Inclination: 

33.94° 

67.172° 

Period: 

140 hrs . 

104 hrs. 

Spin Rate: 

23 R.P.M. 

23.4 R.P.M. 



TABLE II 

COEFFICIENTS OF TELESCOPE RESPONSE FUNCTIONS 


Experiment 

a 

b 

f 

£ 

IMP-III 

-3.104 

0.776 

0.0679 

0.862 

IMP -IV MED, 

o.o ■ 

0.647 

0 . 0434 

0.749 

IMP- IV LED 

-0.272 

0.341 

-0.0945 

0.0645 
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TABLE III 

ANALYSIS PARAMETERS 


Bin No . 

^min 

E max 

E 

4E 

G 

IMP-III 

and IMP-IV MED 

Proton: 




1 

18.7 

29.2 

23.9 

10.5 

3.26 

2 

29.2 

39.7 

34.4 

10.5 

3. 14 

3 

39.7 

50.2 

44.9 

10.5 

2.99 

4 

50.2 

60.6 

55.4 

10.4 

2.84 

5 

60.6 

71.1 

65.8 

10.5 

2.67 

6 

71.1 

81.7 

76 .4 

10.6 

2.52 


IMP-III and IMP-IV MED Alpha Particle: 


1 

18.7 

40.0 

29.4 

21.3 

3.20 

2 

40.0 

60.0 

50.0 

20.0 

2.92 

3 

60.0 

81.5 

70.8 

21.5 

2.62 

IMP-IV 

1 

LED Proton: 

4.2 

6.2 

5.2 

2.0 

.764 

2 

6.2 

9.5 

7.8 

3.3 

.756 

3 

9.5 

13.5 

11.5 

4.0 

.740 

4 

13.5 

19.1 

16.3 

5.6 

.594 

IMP-IV LED Alpha Particle 
1 4.8 

7.0 

5.9 

2.2 

.764 

2 

7.0 

10.0 

8.5 

3.0 

.756 

3 

10.0 

13.7 

11.8 

3.7 

.740 

4 

13.7 

18.9 

16.3 

5.2 

.594 

Note: 

Energies in units 

of MeV/nucleon. 

Geometry factor in units 

of cm2 



TABLE IV 


SEMI-ANNUAL COUNT OF 29 MEV 
PROTON FLUX INCREASES 



Dates 


“3* 

> 10 J 

> 10” 2 

IV 
f— 1 

o 

1 

f -* 1 

IV 

1 — 1 
o 
o 

> 10 1 

Jan. 

- Jun. 

'65 

0 





Jul. 

- Dec. 

'65 

1 

1 




Jan. 

- Jun. 

'66 

4 

1 




Jul. 

- Dec . 

'66 

8 

4 

3 

1 

1 

Jan. 

- Jun. 

'67 

(18)** 

(14) 

(8) 

(5) 

(2) 

Jul. 

- Dec . 

'67 

7 

2 

2 



Jan . 

- Jun . 

'68 

6 

5 

1 

1 


Jul . 

- Dec . 

'68 

(16) 

(12) 

(4) 




TOTALS 


60 

39 

18 

7 

3 


* Units are number/s cm^ sr MeV. 

'* The numbers in parenthesis are normalized to 6 months where data 


was missing. 



Particles / M z -sec-sterad. MeV/Nucleon) 
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Figure 1. Primary differential energy/nucleon spectra of cosmic-ray 
protons and helium nuclei observed near earth near solar 
minimum in 1965. Solid and open symbols represent 
measurements of protons and helium nuclei, respectively. 

1. IMP-3, June-Dee. 1965, Fan et al. (1965a) 

2. IMP-3 and balloon, June 1965, Balasubrahmanyan, Hagge, 
et al. (1966a, b) 

3. balloon, June 1964, Waddington and Freier (1966) 

4. balloon, June 1963, Freier and Waddington (1965) 

5. balloon, June 1964, 1965, Ormes and Webber (1966) 

6. balloon, 1954, McDonald (1958) 

7. balloon, May 1965, Hofmann and Winckler (1966) 

8. balloon, April 1963, Anand et al. (1966) 

Only data representative of solar minimum have been used, 
(from Gloeckler and Jokipii, 1967). 
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Figure 2. Comparison of Deep River neutron counting rate with smoothed 
sunspot number for the current cycle. Mean of previous 19 
cycles superimposed. 
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Figure 4. AE versus E - AE response curves for hydrogen and helium 
isotopes for the scintillator telescope used on IMP-III 
and IMP- IV. 





Figure 5. AE versus E - AE response curves for hydrogen and helium 
for the solid state particle telescope on IMP- IV. 






Figure 6. Schematic diagram of geometry of a generalized circular 

cylindrical telescope of height h and base radii of R x and R £ . 



Figure 7. Schematic of IMP scintillator telescope assembly. 



m„ = 168.4 G 


Figure 8. Schematic of IMP scintillator telescope showing critical 
dimensions and composition of layers. 
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PENETRATION OF B DETECTOR (CM) 

Figure 9. Geometry factor G versus penetration of particle into B 
detector for IMP scintillator telescope. 



17.526 CM. 



END CAP (AL ALY) 
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LIGHT SHIELD (TITANIUM) 
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INSERT (PILOT B PLASTIC 
SCINTILLATOR) 

DETECTOR HOUSING (AL ALY) 

WASHER (SILICONE RUBBER). 

( 2 REQUIRED ) 

300 SQ. MM. RING MOUNTED 
SURFACE BARRIER DETECTOR 
(I MM THICK, 2 REQUIRED) 
WASHER (COPPER, 5 REQUIRED) 
END PLUG (PILOT B PLASTIC 
SCINTILLATOR) 


PHOTOMULTIPLIER TUBE 
( RCA #C-7I5 1 ) 
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TUBE HOUSING (AL ALY) 


GEOMEYRY FACTOR 5 CM -SR 


FILTER CIRCUITRY (PC BOARD) 


POWER SUPPLY, PS -15/ 1 


Figure 10. Schematic of IMP- IV solid state detector telescope assembly. 



0.000229 CM. Ti FOIL 



.0146 CM. Si (A CRYSTAL) 

.00004 G CM” 2 Au (TOP) 

,00004 G CM" 2 AI (BOTTOM) 

R= 0.9525 CM 

.00004 G CM -2 Au (EACH SURFACE) 

0.1000 CM. Si (B| CRYSTAL) 

0.0984 CM Si (B CRYSTAL) 

.00004 G CM" 2 AI EACH SURFACE 


p Tl =4.540 G CM 
p An =19.320 G CM” 3 
/> A| = 2.702 G CM' 3 
p si = 2.330 G CM' 3 


Figure 11. Schematic of IMP- IV solid state detector telescope showing 
critical dimensions and composition of layers. 
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Figure 14. Block diagram of IMP- IV cosmic ray experiment electronics. 
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Figure 15. IMP- III AE versus E-AE matrix showing proton line and 
six energy bins: (1) 18.7 - 29.2, (2) 29.2 - 39.7, 

(3) 39.7 - 50.2, (4) 50.2 - 60.6, (5) 60.6 - 71.1, 
and (6) 71.1 - 81.7 MeV. 
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Figure 17(a). 24 hour average proton fluxes. 
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Figure 17(b). 24 hour average proton fluxes. 
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Figure 17(c). 24 hour average proton fluxes. 


Ill 



FLUX NO. / SEC CM‘ SR M«V 


10 


10 ' 


2'i 


10 


10 


.-2 


!’*'’•**** C j 1 ! i 

! ! . ! ! I I! 

1 ! ! i i i it 

• 

i t 
i i 
i i 
i t 
i i 
i t 
i i 

i i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 

i r 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 
1 1 
1 1 
1 1 

1 1 
1 1 
1 1 

| 

i i 
\ i 
t i 
i i 
i i 
* i 
\ i 
i i 
j i 

1 

1 

1 

t 

1 

1 

1 

! i 

1 1 

1 1 

i : 

1 1 
f 1 
1 1 

1 1 
1 1 
1 1 
1 1 

• 

i j 

i i 

i i 

i i 

i i 

i i 

i i 

1 

1 

1 

1 

1 

1 

1 

1 

" 1 1 

1 1 

1 1 

1 | X 

• t x| 

» 1 1 X 

1 0 0 

1 1 

1 1 
1 1 
1 1 
1 1 

1 1 
1 1 
1 1 
1 1 

i C |0 || » 

*t i m i x *i x 

XC |C X| lo 

It 1 1 1 o xo |0 1 0 x| 



0 1 1 IC Q || 1 la II 

X 1 X| |C X) 

XO XO xo X| K X 1 1 Id 

I x| | 

i i 

1 c 0 0 O MO 

II 1 1 • II 

1 '1 ' 

1 1 1 1 O 1 * 1 I . 

! i 

lo 1 
1 1 


— — — XO— < 

I 

xc 
x| I 
lo 


■„ — - 4 • — — < 


-*0 1 


1 1 * = 12 MeV 

1 1 
1 1 

' 1 x * 29 MeV 

1 t 
»••• 

! ! o = 60 MeV 

i i 

i i 

i i 

i i i i 

i. i i i 

i . \ i i 

1 t 
1 1 
t 1 
1 « 
1 1 
1 I 
1 1 
1 1 

i i i ii 

i i i i r 

* ‘ * i i it 

! i ii 

! 1 1 1 * 

! ! 1 1 i i 

! 1 i i 

i i i i ii 

t i i i 

1 x t I 

1 X 1 .1 1 

1 X | 1 1 

1 1 X t 1 

1 a X| 0 | | 

o x X I | 1 

1 1 
1 1 
1 1 
1 1 
1 • 
1 xt 
1 t 

1 C | X 1 1 

1 O 1 OX, 1 

»t ! 

1 1 0 1 1 0 |o xo XO 1 } 

1 t 0 0 X | II 0 0 

1 1 III 

i i . i 

XX II 

1 It 0 t 

O’ o 1 1 

t 1 t 


JUN 1966 


JULY 1966 


K> 


-2 

10 


K> 


i i r i r i 

i : : i i . i 

— 

t 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

t 

i 

1 

t 

1 

1 

1 

1 

1 

1 

1 

t 

1 1 t 

1 1 t 

1 1 1 

1 1 t 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

t 

1 

1 

1 

1 

1 

1 

J 

1 

1 

1 

1 

t 

1 

1 

1 

t 

1 

1 

1 

1 

1 * 
t 
1 
1 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

t 

1 

» 

1 

1 

1 

t 

1 

1 

t 

1 

t 

i 

i 

i 

i 

i 

j 

i 

i 

i 

i 

i 

i 

i 

I 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

t 

t 

1 

1 

1 

t 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 1 1 

1 t » 

1 1 1 

1 1 1 

1 1 I 

1 t t 

1 I I 

1 1 1 

1 

t 

1 

1 

« 

1 

J 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

t 

1 

1 

1 

t 

1 

1 

t 

1 

t 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I X X 

1 
1 
t 

t 0 

j, 

i 

i 

i 

t 

] 

i 

i 

i 

* 

i 

i 

x| 

1 

x 1 

0 1 X 

0 

1 X 

1 

t 

1 

i 

1 

1 

1 

1 

1 

1 
1 
1 
1 
1 
1 
1 
X J 

x 1 

1 

1 

1 

1 

1 

1 

t 

1 

1 

1 I t 

1 1 1 

1 1 1 

t 1 1 

f t • 

1 1 t 

1 1 1 

1 1 1 

t 

1 

1 

1 

i 

1 

1 

1 

t 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 o 

1 

1 

x | 
1 
1 
1 

i 

i 

1 X 

i 

1 a 

C } 

1 1 > 

1 C 

1 X 

1 

1 

I 

1 

1 

1 

1 X 

c j 

0 x 

1 0 

j 

1 

1 

1 X 

x 1 

t 1 1 

1 1 t 

1 x It 

1 t 1 

1 X I • 

1 1 
1 > 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

- 1 

1 

1 

1 1 

1 

1 

t 

1 

1 c 
1 

1 c 

1 

1 0 

xi 

j X X 

1 

x | o 

! o 

l 

1 

1 

0 1 

1 loo || 

x 1 X II 

x i ii 

1 tit 

I 1 C*| »| 

1 c 1 e |c 
1 1 1 

IC to 

X* *1 X 
1 1 

1 0 

x ( x| c x| o 

1 xo 1 X O X 1 

C II J | 

1 

x| 0 | 

X 1 | xo 

II 0 1 j 

0 t X 

1 XI 1 XO 
xo lo XO || X| 
1 H II lo 

1 I* 

X It 1 

|0 | II | 0 xo 
1 *| *0 X| 1 
O II 1 

1 

1 

t 

1 

1 

1 

1 

1 

1 

1 

1 1 O 

1 c 
1 

1 

o I 

1 
1 
1 

1 

1 

1 

1 

o o 1 xl 1 

1 O xo || 

1 II o a l 

!! i j i 

0 

5 

to 

15 

20 

25 3 

0 

0 

5 

to 

15 

20 

25 30 


AUG 1966 


SEPT 1966 


Figure 17(d). 24 hour average proton fluxes. 
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Figure 17(f). 24 hour average proton fluxes. 
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Figure 17 (i). 24 hour average proton fluxes. 
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Figure 17(j). 24 hour average proton fluxes. 
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Figure 18. IMP-III and IMP-IV monthly averages of proton fluxes and 
monthly averages of the Deep River hourly neutron rate. 
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Figure 21. Quiet time alpha particle energy spectra. 
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Figure 24. Time history of the 24 hour relative particle fluxes from 
the IMP-IV discriminator level 6. 
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Figure 25. Representative 96 hour proton and alpha particle energy 
spectra with the least squares fit to the two component 
function. 
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Figure 26 (b) , Averages of alpha particle-to-proton ratios at 
discrete energies from 24 May 1967 to 20 August 
1968 showing standard deviations. 
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Figure 27. 76.4 MeV versus 5.2 MeV proton fluxes. 
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Figure 28. 5.2 MeV proton flux versus the energy at which minimum 
flux occurs. 
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Figure 29. 76.4 MeV proton flux versus the energy at which minimum 
flux occurs. 
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Figure 31. The computed power law indices versus the S.2 MeV proton 
flux. 
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Figure 33. 74.2 MeV/nucleon versus 5.9 MeV/nucleon alpha particle 
fluxes. 
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Figure 34. 5.9 MeV/nucleon alpha particle flux versus the energy/ 
nucleon at which the minimum flux occurs. 
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Figure 35. 74.2 MeV/nucleon alpha particle flux versus the energy/ 
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Figure 36. The minimum alpha particle flux versus the energy/nucleon 
at minimum flux. 
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Figure 37. 


The computed power law indices versus the 5.9 MeV/nucleon 
alpha particle fluxes. 
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Figure 38. The computed power law indices versus the 74.2 MeV/nucleon 
alpha particle fluxes. 
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Figure 39. The 24 hour average proton fluxes at 5.2 MeV for IMP-IV. 

The times large confirmed solar flares occrrred are marked 
F. Central meridian passage of calcium plage regions are 
marked by the plage region number. Recurring regions are 
connected by horizontal lines. 
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Figure 40. 10% variable time averages of proton fluxes for 3 months 
of imp- IV data showing solar proton increases caused by 

corotating regions. 
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Figure 41. 6 hour integral proton fluxes for energies >20 MeV from 
IMP-III. The times of meridian transit for returning 
calcium plage regions are marked with the region numbers 
Flare times are marked with an F. 








